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1.  Introduction 


The  purpose  of  this  research  is  to  develop  a  multivariate  autore¬ 
gressive  (MVAR)  climatological  model  to  be  used  for  short-term  forecasting 
(3-12  hours)  of  various  atmospheric  variables  over  a  limited  area  in  a 
tactical  situation.  The  atmospheric  variables  to  be  forecast  are  temperature, 
visibility,  ceiling  height,  height  of  the  first  cloud  layer,  and  the  u-  and 
v-components  of  the  wind.  The  tactical  situation  facing  the  forecaster  is 
this:  all  of  his  communications  are  cut  off  and  he  must,  using  only  a  smal’ 
computer  with  limited  storage  area,  make  3-hourly  forecasts  of  the  afore¬ 
mentioned  variables  out  to  12  hours  over  an  area  on  the  order  of  100  km 
square.  Using  an  MVAR  model  the  forecaster  can  not  only  make  the  necessary 
forecasts,  but  confidence  intervals  about  the  forecast  values  can  also  be 
computed  to  aid  in  any  decision-making  processes  based  on  these  forecasts. 

In  the  next  section  a  theoretical  description  of  an  MVAR  forecast 
model  is  presented.  In  such  a  model  the  forecast  value  of  a  variable 
(predictand)  is  assumed  to  be  a  function  of  present  and  past  observations  of 
that  variable  as  well  as  other  predictor  variables.  The  relationships 
between  the  predictand  and  the  predictors  are  carried  within  the  coefficient 
matrices  of  the  model  which  are  determined  from  the  past  history  of  observations. 
Embedded  within  these  coefficient  matrices  will  be  the  effects  of  complex 
terrain  upon  the  inter-relationships  among  the  predictand  and  predictors. 

The  data  used  in  this  study  consisted  of  five  years  of  3-hourly 
observations  of  the  six  variables  to  be  predicted  as  well  as  one  predictor  variable, 
the  dewpoint  temperature,  for  the  following  north  German  stations:  Hannover, 
Bremen,  Braunschweig ,  Boizenburg,  Magdeburg,  Weissen,  and  Wernigerode.  In 
Section  3  the  preliminary  data  analysis  is  described  in  which  the  distribution 
and  various  statistical  properties  of  the  variables  were  determined.  From 
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the  re-  t'-  of  this  analysis  it  was  decided  that  the  c'oud  he i <;h * 
ibility  variables  be  transformed  ■>>  such  a  way  that  the  variance  pr< 
during  low  visibility  and  ceiling  conditions  be  emphasized  so  tor* 
conditions  could  be  better  forecast. 

Finally,  in  Section  4  various  MVAR  forecast  models  wr  re¬ 
determined  using  the  first  three  years  of  data  (  1972-1974'  .  W”-/-  • 

to  the  limited  data  storage  available  to  the  forecaster  it  wa  .  a> 
the  final  models  could  only  util  ize  values  of  the  predictor-,  t*  ■ 
period  preceding  the  forecast  time.  It  was  found  that  under  th>  t.: 
constraints  the  best  forecasts  could  be  obtained  using  six  7-parameter 
MVAR  models.  Each  model  produces  a  forecast  for  a  particular  variat'-’  - 
the  observations  at  the  seven  stations  as  parameters.  The  variab’o.  *  ■ 
one  can  theoretically  expect  to  forecast  best  are  temperature  arc  fr" 
and  v-components  of  the  wind  with  about  Q5,  75  and  60  percent  o*  the 
variance,  respectively,  explained  by  the  model.  The  percent  variance  cni 
expect  to  explain  ranges  from  about  45  to  70  for  visibility  and  from  -i‘: 

30  to  60  for  the  cloud  height  variables.  Using  the  data  from  1976  tk,e  vi ' 
models  were  tested  and  the  error  statistics  from  these  actual  forecasts  w<  * 
found  to  agree  quite  well  with  theory. 


! 
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2.  Multivariate  Autoregressive  Forecast  Model 

In  this  section  the  theory  behind  a  multivariate  autoregressive 
(MVAR)  forecast  model  will  be  outlined.  Suppose  that  one  has  collected  m 
time  series  consisting  of  n  observations  each  for  m  different  variables. 
Further  assume  that  the  value  of  the  sample  mean  has  been  subtracted  from 
each  of  these  observations.  Using  these  observations  one  would  like  to 
develop  a  forecast  model  such  that  future  values  for  m,  <  m  of  these 

i  — 

variables  may  be  predicted,  given  the  present  and  a  certain  number  of  past 
values  of  these  variables.  For  a  particular  time  i  the  observations  tor  the 
m  variables  are  denoted  by  the  m-dimensional  column  vector,  ,  where  the 
first  elements  of  X.  belong  to  the  time  series  to  be  predicted  while 
the  last  m  -  mj  belong  to  those  series  to  be  used  to  aid  in  the  forecast. 

The  p-th  order  MVAR  model  is: 

♦  AjXf.j  +  ...  ♦  Vh-p  =  zi  (1) 

where  the  A's  aremx  m coefficient  matrices  and  is  an  m-dimensional  white 
noise  column  vector.  Such  a  model  would  use  the  present  observations  (X-j_j) 
and  p-1  past  observations  in  order  to  predict  the  values  of  the  variables 
one  interval  in  the  future  (X^).  The  variance  of  the  white  noise  process 
(Z-j)  represents  the  one-step  prediction  error  variance  of  forecasts  made 
with  this  model.  Two  things  must  be  determined, using  the  collection  of  n 
observation  vectors,  before  an  MVAR  forecast  model  can  be  developed.  First, 
the  proper  order  model  must  be  selected  and  then  the  coefficient  matrices 
for  that  order  model  must  be  computed. 

The  procedure  outlined  in  this  section  is  the  multivariate  general¬ 
ization  given  by  Whittle  (1963)  of  the  recursive  method  developed  by  Durbin 
(1960)  for  the  fitting  of  univariate  autoregressive  models  of  successively 
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increasing  order.  Except  for  the  inclusion  of  the  Akaike  FPE  criterion 
(Akaike,  1971),  it  is  identical  to  that  presented  by  Jones  (1964).  The 
Akaike  FPE  parameter  is  an  estimator  of  the  one-step  prediction  error  of 
the  MVAR  process.  The  use  of  the  FPE  criterion  permits  one  to  find  the 
order  model  with  the  smallest  one-step  prediction  error.  The  analysis 
procedure  can  be  described  simply:  first,  L  +  1  MVAR  models  whose  order 
successively  increase  from  zero  to  L  are  fitted  to  the  n  m-dimensior  ’  ’ 
observation  vectors  using  the  recursive  method  to  be  detailed  below, 
for  each  order  model  a  value  of  the  Akaike  FPE  parameter  is  computed,  an., 
finally,  an  MVAR  model  whose  order  is  that  for  which  the  minimum  FPE  w’S 
found  is  fitted  to  the  data  using  the  same  recursive  method.  This  ‘is  the 
model  that  one  would  use  for  prediction. 


The  first  step  in  the  analysis  method  is  the  calculation  of  the 


lag  sums 


G  =  .  L 
P  l=p+l 


i=p+l  XiXi-P’  P=0,1’2’"-’L’ 


6-P  =  GP 


where  denotes  the  transpose  of  the  m  x  m  matrix  Gp.  In  the  following 
equations  the  p-th  order  residual  matrices  for  the  forward  and  backward 
autoregressions  are  denoted  by  Sp  and  Sp,  respectively.  The  k-th  coeffi¬ 
cient  matrices  for  the  p-th  order  forward  and  backward  autoregressions  are 
denoted  by  AP  and  respectively.  The  determinant  of  the  m^  x  m^  sub¬ 
matrix  in  the  upper  left-hand  corner  of  Sp  is  denoted  by  (Sp,mj|. 


Initial ization: 


a;  =  -  g_i 


FPEo 


in+lviH-. 

(n^T}  1 


i SQ ,m  j j  . 


For  p  =  1,2,  . . . ,  L-l , 


rina1 ly. 


sp  =  Go  *  fl?  G-i  *  *  Ap  G-P 


S  =  G  +  A?  G.  +  ...  +  A?  G 
Poll  pp 


,n+pm+l  ,m. 


FPEP  *  <ST5STT>  1  iVll 

Cl  ■  -<vi  * A?  v-tA!6i>  v1 

^  6-P  *  ---  *  *s  «-i> 


Ad+1  =  Ap  +  AP+1  ^ 

k  Hk  Mp+1  p+i-k 

AP+1  =  aJ  +  aP+1  ap 

k  k  p+i  p+l-k 


SL  -  G0  +  AE  G_j  *  ...  *  AlL  G_l 


k  =  1,2,  ...»  p. 


FPE, 


_  /  n+Lm+1  \m-i  i<-  : 

l  ln=rsrr  'W- 


Thus,  it  one  had  found  that  FPE^  had  been  a  minimum  and  had  fitted  at  k-th 
order  MVAR  model  to  the  data,  the  prediction  for  given  X^-,  ...,  X^ 
would  be 


h  =_  .  xi-i 


*2  xi-2 


-  A,  X.  , 
k  i-k 


(2) 


o 


At  this  point  the  values  of  the  sample  means  for  the  variables  to  be  predicted 


would  be  added  back  to  the  X-  vector  to  give  the  final  prediction  value  fo*- 
each  variable. 

Finally,  one  can  determine  the  quality  of  the  predictions  from  such 
a  model  by  computing  the  prediction  error  covariance  matrices.  For  a  k-t'n 
order  model  the  one-step  prediction  error  covariance  matrix  is 

w1  =  1  c 

Vk  n-mk  * 

Successive  predictions  can  be  made  using  (2)  by  merely  reolacinq  ebserv-  ■ 
values  by  predictions  as  one  steps  further  into  the  future.  The  fo'low'rc 
recursion  is  used  to  find  the  i-steD  prediction  error  covariance  "iatv"'^, 

VJ,  when  using  repeated  predictions:  1 

5,  •  -  Ai 

Bj  *  -  (Aj  +  +  •••  *  3>i  V 

and 

Vj  -  V4-1  4-  Bj.i  V1  Bj.j  . 

Once  this  matrix  has  been  obtained,  its  main  diagonal  consists  of  the  error 
variances  for  the  variables  to  be  predicted.  These  can  then  be  used  to 
determine  confidence  intervals  to  be  placed  about  the  forecast  values. 


3.  Preliminary  Data  Analysis 


Before  an  attempt  was  made  to  determine  any  MVAR  models,  a  pre¬ 
liminary  data  analysis  was  performed  using  the  3-hourly  observations 
(00Z,  03Z,  etc.)  collected  during  1972- 1975  for  the  fo’lowing  stations- 
Hannover,  Bremen,  Braunschweig,  Boizenburg,  Magdeburg,  Weissen,  and 
Wernigerode.  First,  the  distributions  of  the  six  variable"  c  be  oredicted 
(temperature,  u-  and  v-compcnents  of  the  wind,  visibility,  ceiling  height, 
and  height  of  the  first  cloud  layer)  were  determined.  In  Figure  i  the 
distribution  of  the  1972-1975  temperature  and  u-wind  observations  fo1" 
Wernigerode  are  displayed  with  their  sample  means  denoted  by  a  star.  The 
distributions  shown  here  are  typical  cf  the  temperature  and  u-  and  v-winc 
observations  for  all  seven  stations  used  in  this  study.  As  can  be  seen  in 
Figure  1,  these  variables  appear  to  be  quite  normally  distributed. 

On  the  other  hand  the  distributions  found  for  visibility  and  the 
cloud  height  variables  were  far  from  normal.  In  Figure  2a  the  distribution 
found  for  the  visibilities  observed  at  Wernigerode  is  shown.  This  distri¬ 
bution,  which  is  typical  of  those  found  for  the  visibility  and  cloud  height 
variables  at  all  seven  stations,  is  roughly  rectangular  with  a  sample  mean 
of  just  over  9  km.  However,  since  an  MVAR  model  is  designed  to  predict 
deviations  of  the  variables  about  their  sample  means  and  since  for  these 
three  variables  the  sample  means  are  much  larger  than  the  low  visibility 
and  low  ceiling  situations  that  one  would  like  to  be  able  to  predict,  the 
following  transformations  were  performed  upon  these  variables  in  order  to 
emphasize  the  variance  of  the  low  visibility  (ceiling)  situations: 


V 2  =  exp(-Vl/2000) 


NO.  (Jh  OBSERVATIONS 


8000 


6000  \- 


VJSieiL! FT  (KM) 


Figure  2a.  Distribution  of  1972-1975  visibility  observations  for 
Werni gerode. 
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C2  =  exp(-Cl/1000) ,  (4) 

where  VI  and  Cl  are  the  observed  visibility  and  cloud  height,  respectively, 
in  meters.  In  Figure  2b  we  can  see  the  effect  of  this  transformation  upon 
the  visibility  distribution  for  Wernigerode.  The  sample  mean  of  the  trans¬ 
formed  variable  now  represents  a  visibility  of  only  2.4  km.  Furthermore, 
an  observation  of  visibility  less  than  1  km  will  contribute  more  to  the 
variance  than  an  observation  of  unlimited  visibility  since  its  deviation 
from  the  sample  mean  will  be  larger.  By  expanding  the  scale  for  low  vi  a¬ 
bilities  and  decreasing  the  scale  for  high  visibilities,  this  trans forma tie 
permits  more  precise  forecasts  of  the  poor  visibility  situations. 

Next,  monthly  and  hourly  means  were  computed  for  all  variables  at 
all  stations  using  the  1972-1975  observations.  Table  1  summarizes  the 

Table  1.  Variance  explained  by  monthly  and  hourly  means  for  Hannover 
_ (1972-1975). _ 

Total  Variance  Explained  by: 

Variable  Variance  Monthly  Means  Percent  Hourly  Means  Percent 


First  Cloud  laye 

r  Ht.*  .0791 

.0052 

6.6 

.0036 

4.6 

Ceiling  Ht.* 

.0843 

.0070 

8.3 

.0028 

3.3 

Temperature 

53.49 

32.13 

60.1 

4.70 

8.8 

Visibil ity* 

.0331 

.0021 

6.3 

.0013 

3.9 

u-wind 

13.63 

.79 

5.8 

.062 

0.5 

v-wind 

6.08 

.52 

16.4 

.32 

0.5 

*  Indicates  transformed  variable 

• 

results  of  these 

computations  for 

Hannover, 

which  again 

are  typical 

of 

those  found  for  the  other  stations.  We  can  see  that  with  the  exception 
of  monthly  averages  for  temperature,  only  a  small  percent  of  the  total 
variance  for  the  six  variables  is  explained  by  the  annual  and  diurnal  cycles 
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Figure  2b.  Distribution  of  transformed  1972-1975  visibility  observations 
for  Wernigerod. 
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Typically,  about  60  percent  of  the  temperature  variance  i^.  explained  by 
the  annual  cycle.  Plots  of  the  hourly  and  monthly  means  for  Hannover  are 
displayed  in  Figure  3.  We  can  see  that  the  amplitude  of  the  annual  wave  is 
about  8°C  while  that  of  the  diurnal  wave  is  about  3°C.  Since  all  seven 
stations  display  this  pronounced  annual  wave  for  temperature  and  it  explains 
a  significant  amount  of  the  variance,  its  effect  will  be  removed  from  the 
temperature  data  along  with  that  of  the  diurnal  cycle  and  the  sample  mean 
before  any  MVAR  models  are  determined  for  that  variable.  One  must  be  care¬ 
ful  when  attempting  to  fit  an  MVAR  model  to  data  which  are  highly  correlated 
since  unstable  processes  can  be  produced.  In  any  case,  it  is  the  deviations 
about  these  very  regular  cycles  that  we  are  interested  in  forecasting,  and, 
thus  it  is  these  deviations  which  we  will  attempt  to  model. 
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fEMPERBTURE 


HANNOVER  (1972-1975) 
O  -  MONTHLY  MEANS 

A  -  HOURlT  means 


UOZ  03  Z  U6Z  09Z  \21  i5Z  i  6 Z  Li  2 


Figure  3.  Monthly  and  hourly  means  computed  from  1972-1975  Hannover 
temperatures. 
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4.  Application  and  Results 


Using  observations  collected  during  1972-1974  the  MVAR  modeling 
process  described  in  Section  2  was  applied  in  several  different  ways  in 
order  to  determine  the  best  forecasting  technique.  Given  the  tactical 
constraint  of  limited  computer  size  and  storage,  we  will  only  consider 
MVAR  forecast  models  of  order  nine  or  less.  Thus,  only  a  one  day  history 
of  observations  of  the  predictor  variables  need  be  stored  at  a  time.  This 
tactical  consideration  also  limits  the  number  of  variables  to  be  used  by 
the  MVAR  models  (not  only  those  to  be  predicted  but  also  those  to  aid  in 
the  prediction).  For  example,  if  an  MVAR  model  were  developed  which  used 
all  of  the  observed  variables  for  all  stations,  then  it  would  possess  49 
variables.  If  the  model  were  of  order  nine,  then  nine  49  x  49  coefficient 
matrices  would  have  to  be  stored  and  used  in  the  forecast  computations  a1  one 
with  nine  49  x  1  observational  vectors  from  the  past  24  hours.  This  would 
require  almost  100K  bytes  of  storage  and  a  comparably  large  number  of 
computations  required  to  make  the  forecasts.  On  the  other  hand,  a  model  of 
order  nine  with  only  seven  variables  would  require  about  50  times  less  space 
and  computation  time. 

The  first  type  of  MVAR  model  to  be  tested  used  the  six  variables 
to  be  predicted  (temperature,  u-  and  v-wind  components,  visibility,  ceiling 
height,  and  height  of  the  first  cloud  layer)  along  with  the  dew-point  tem¬ 
perature  for  one  station  at  a  tine.  The  analysis  procedure  outlined  in 
Section  2  was  applied  using  the  1972-1974  observations  to  calculate  the 
lag  sum  matrices.  The  model  order  was  restricted  to  be  no  greater  than 
nine  and  this  maximum  was  chosen  in  every  case  except  for  Braunschweig  where 
an  eighth  order  model  possessed  the  smallest  Akaike  FPE  parameter.  The 


IB 


7x7  coefficient  matrices  for  each  resulting  station  model  were  determined 
along  with  estimates  of  the  one-step  prediction  error  variance  for  the  six 
variables  to  be  forecast.  We  will  use  these  estimates  in  order  to  determine 
the  quality  of  the  various  MVAR  forecast  models  tested  here.  Table  2  summar 
izes  the  results  for  the  7  station  models.  We  can  see  that  using  such  a 


Table  2.  Estimates  of  one-step  prediction  error  variance  and  percent 

variance  explained  (in  parenthesis)  for  the  station  MVAR  models 
_ determined  from  1972-1974  observations. _ 

Height  of 

Model  1st  Cloud  Ceiling 


Station 

Order 

Layer* 

Height* 

Temperature 

Visibility* 

u-wind 

v-wind 

Hannover 

9 

.0374 

.0436 

2.55 

.0164 

2.96 

2.40 

(54.7) 

(48.6) 

(94.8) 

(52.2) 

(78.8) 

(63.0) 

Bremen 

9 

.0367 

.0436 

2.93 

.0194 

3.31 

2.99 

(52.8) 

(45.5) 

(93.8) 

(42.4) 

(78.4) 

(67.3) 

Boizenburg 

9 

.0181 

.0308 

3.75 

.0252 

4.78 

3.47 

(32.2) 

(35.0) 

(92.4) 

(46.8) 

(63.9) 

(48.4) 

Braunschweig 

8 

.0296 

.0333 

3.12 

.0112 

2.82 

2.46 

(53.2) 

(49.8) 

(92.7) 

(55.0) 

(75.9) 

(59.1) 

Magdeburg 

9 

.0260 

.0314 

3.58 

.1400 

3.01 

3.17 

(36.3) 

(38.8) 

(93.3) 

(66.7) 

(72.1) 

(52.0) 

Wernigerode 

9 

.0223 

.0264 

3.9 

.0332 

3.97 

4.16 

(38.4) 

(42.5) 

(92.1) 

(64.6) 

(67.5) 

(46.0) 

Weissen 

9 

.0360 

.0391 

4.15 

.0186 

5.01 

3.06 

(33.1) 

(38.8) 

(92.4) 

(50.4) 

(67.5) 

(52.9) 

*Indi cates'  transformed  variable. 


model  we  can  best  forecast  temperatures  (from  92.1  to  94.8  percent  variance 
explained)  and  are  least  able  to  forecast  the  height  of  the  first  cloud 
layer  (from  32.2  to  54.7  percent  variance  explained).  One  can  roughly 
expect  to  explain  40%,  55%,  55%,  and  70%  of  the  variance  for  ceiling  height, 
visibility,  v-wind,  and  u-wind,  respectively. 


The  next  type  of  MVAR  model  to  be  investigated  utilized  the  obser¬ 
vations  for  a  particular  variable  at  all  seven  stations.  Again  the  model 
order  was  restricted  to  nine  or  less  and  the  MVAR  analysis  procedure  was 
applied  to  the  1972-1974  observations  for  the  six  different  variables  to 
be  forecast.  In  Table  3  the  estimates  of  the  one-step  prediction  error 
variance  and  the  percent  variance  explained  for  the  6  variable  models  are 
displayed.  In  this  case,  as  discussed  in  Section  3,  the  monthly  and  hourly 


Table  3.  Estimates  of  one-step  prediction  error  variance  and  percent  variance. 

explained  (in  parenthesis)  for  the  variable  MVAR  models  determined 
from  1972-1974  observations. 


Model 

Variable  Order  Hannover  Bremen  Boizenburg  Braunschweig  Magdeburg  Wernl ge-rode  Wejssen 

Height 
of  1st 


Cloud 
Layer  * 

9 

.0359 

(56.5) 

.0369 

(52.5) 

.0183 

(31.5) 

.0261 

(58.7) 

.0244 

(40.2) 

.0224 

(38.1) 

.0339 

(37.0) 

Ceiling 
Height  * 

9 

.0419 

(50.6) 

.0440 

(45.0) 

.0298 

(37.1) 

.0293 

(55.9) 

.0282 

(45.0) 

.0263 

(42.7) 

.0367 

(42.6) 

Tempera¬ 

ture 

9 

2.24 

(95.4) 

2.34 

(95.0) 

2.63 

(94.7) 

2.02 

(95.3) 

2.35 

(95.6) 

2.82 

(94.3) 

3.19 

(94.1) 

Visi¬ 
bility  * 

9 

.0157 

(54.2) 

.0187 

(44.5) 

.0241 

(49.2) 

.0097 

(61.0) 

.0132 

(68.6) 

.0335 

(64.3) 

.0178 

(52.5) 

u-wind 

9 

2.51 

(82.0) 

3.18 

(79.2) 

3.85 

(70.8) 

2.28 

(80.5) 

2.33 

(78.4) 

3.66 

(70.0) 

3.88 

(74.8) 

v-wind 

9 

2.00 

(69.1) 

3.08 

(66.3) 

2.99 

(55.6) 

2.03 

(66.3) 

2.51 

(62.0) 

3.72 

(51.7) 

2.48 

(61.8) 

*  Indicates  transformed  variable. 


means  were  removed  along  with  the  sample  means  for  the  temperature  observations. 
A  comparison  of  Tables  2  and  3  indicates  that  in  almost  every  case  the  percent 
variance  explained  by  the  variable-at-a-time  MVAR  models  is  greater  than  that 
for  the  station-at-a-time  models.  The  greatest  improvement  is  noted  for  the 


v-wind  predictions  where  the  percent  variance  explained  is  increased  by  as 
much  as  10%.  For  the  variable- at- a- time  models  we  can  expect  to  explain 
approximately  95%,  75%,  60%,  55%,  45%,  and  45%  of  the  variance,  respecti vely , 
for  temperature,  u-wind,  v-wind,  visibility,  ceiling  height,  and  height  of 
the  first  cloud  layer. 

For  the  first  two  types  of  MVAR  models  tested  here  the  model  order 
was  limited.  However,  for  both  types,  MVAR  models  were  determined  in  which 
the  maximum  order  permitted  was  30.  In  these  cases  MVAR  models  whose  order 
ranged  from  25  to  30  were  found  to  possess  the  minimum  value  of  the  Akaike 
FPE  parameter.  In  every  case  though  the  reduction  of  the  one-step  prediction 
error  variance  over  that  of  the  ninth  order  models  was  negligible.  Thus, 
the  limitation  of  the  model  size  required  by  the  tactical  situation  has  no 
detrimental  effect  upon  the  quality  of  the  forecast  models  produced. 

A  final  MVAR  model  was  examined  in  which  the  number  of  variables 
was  21,  consisting  of  the  aforementioned  seven  variables  for  the  stations, 
Hannover,  Bremen,  and  Braunschweig.  Again  the  model  order  was  limited  to 
nine  or  less  and  a  ninth  order  model  was  chosen  by  the  analysis  procedure. 
Table  4  summarizes  the  results  for  this  particular  model. 


Table  4.  Estimates  of  one-step  prediction  error  variance  and  percent 
variance  explained  (in  parenthesis)  for  the  three-station 
_ MVAR  model  determined  from  1972-1974  observations. 


Station 

Height  of 
1st  Cloud 
Layer* 

Ceiling 

Height* 

Temperature 

Visibility* 

u-wind 

v-wind 

Hannover 

.0341 

(58.7) 

.0391 

(53.9) 

2.10 

(94.7) 

.0152 

(55.7) 

2.46 

(82.3) 

1.97 

(69.6 

8 remen 

.0348 

(55.2) 

.0416 

(47.9) 

2.31 

(95.1) 

.0182 

(46.0) 

2.97 

(80.6) 

2.87 

(68.6 

Braunschweig 

.0248 

(60.8) 

.0282 

(57.5) 

2.19 

(94.9) 

.0098 

(60.6) 

2.24 

(80.8) 

1.98 

(67.1 

*  Indicates  transformed  variable. 


Comparing  Table  4  with  Table  3  one  can  see  that  this  21-variable  model  is 
only  slightly  better  than  the  three  respective  7-variable  models.  Therefore, 
since  the  variable-at-a-time  models  can  be  run  using  about  one-tenth  the 
computer  space  and  time  and  since  there  is  negligible  improvement  to  be 
gained  from  the  larger  model,  they  have  been  chosen  as  the  best  MVAR  fore¬ 
cast  model  to  be  used  for  short-range  predictions  in  a  tactical  situation. 

We  have  seen  in  this  section  how  well,  based  on  the  MVAR  theory 
described  in  Section  2,  we  can  expect  to  forecast  the  six  meteorological 
variables  of  interest  in  this  study.  Using  the  six  variable-at-a-t'me 
models,  whose  expected  performances  are  outlined  in  Table  3,  and  observa¬ 
tions  collected  during  1976,  a  number  of  MVAR  forecasts  we*-e  made  and  com¬ 
pared  with  the  observations  valid  at  the  forecast  time.  Assuming  that  the 
MVAR  forecast  models  are  unbiased,  theoretical  estimates  of  the  root  mean 
square  errors  (RMSE's)  for  the  one-step  through  four-step  predictions  are 
obtained  by  simply  taking  the  square  root  of  the  one-step  through  four-step 
prediction  error  variances.  These  theoretical  RMSE's  are  then  compared  with 
the  actual  RMSE's  computed  from  3-,  6-,  9-,  and  12-hour  forecasts  made  using 
the  six  MVAR  models  (developed  using  1972-1974  data)  upon  1976  observations. 

Tables  5  and  6  display  the  results  of  48  MVAR  forecasts  using  the 
ninth  order  variable-at-a-time  models  for  the  transformed  variables  height 
of  the  first  cloud  layer  and  ceiling  height,  respectively.  We  can  see  in 
Table  5  that  in  virtually  every  case  the  RMSE's  computed  from  the  actual 
MVAR  forecasts  were  smaller  than  those  expected  from  theory  and  in  no  case 
were  they  larger.  The  RMSE's  determined  from  the  ceiling  height  forecasts 
shown  in  Table  6  agree  quite  closely  with  their  theoretical  counterparts 
for  four  stations  (Hannover,  Bremen,  Braunschweig,  and  Weissen)  and  are 


Table  5.  Comparison  of  theoretical  forecast  RMSE's  (T)  with  those  computed 
from  actual  MVAR  forecasts  (C)  using  197?  transformed  height  of 
_ first  cloud  layer  data.  _  _ 


RMSE  ’  s 


No 

.  of 

3- 

■hour  FCST 

6: 

■hour  FCST 

9: 

■hour  FCST 

12-hour  I 

:CST 

Station  Forecasts 

C_ 

T 

C_ 

_T_ 

r 

T _ 

c__ 

T  ' 

Hannover 

48 

0. 

12 

0. 

19 

0. 

.18 

0.23 

0. 

.19 

0. 

25 

0. 

20 

0, 

.26 

Bremen 

48 

0. 

11 

0. 

19 

0. 

.16 

0.23 

0. 

.17 

0. 

25 

0. 

20 

0. 

.26 

Boizenburg 

48 

0. 

10 

0. 

14 

0. 

.11 

0.15 

0. 

.12 

0. 

15 

n 

12 

0. 

.15 

Braunschweig 

48 

0. 

13 

0. 

16 

0. 

,14 

C.  19 

0. 

,20 

0. 

21 

C. 

17 

0, 

.22 

Madgeburg 

48 

0. 

09 

0. 

16 

0. 

15 

0.17 

0. 

.11 

0. 

18 

0. 

14 

r 

■ 

.  19 

Wernigerode 

48 

0. 

08 

0. 

15 

0. 

.10 

0.16 

0. 

.10 

0. 

17 

0. 

09 

0. 

.18 

Weissen 

48 

0. 

17 

0. 

18 

0. 

20 

0.20 

0. 

.21 

0, 

21 

'■  j  . 

22 

0, 

.22 

Table  6.  Comparison  of  theoretical  forecast  RMSE's  (T)  with  those  computed 
from  actual  MVAR  forecasts  (C  )  using  1976  transformed  ceiling 


height  data. 


RMSE's 


No, 

.  of 

3-hour 

FCST 

6-hour 

FCST 

9-hour 

FCST 

12-hour 

rCS" 

Station  Forecasts 

C 

T 

r 

T 

r 

T 

r 

T 

Hannover 

48 

0.20 

0.20 

0.24 

0.24 

0.28 

0.25 

0.25 

0.27 

Bremen 

48 

0.22 

0.21 

0.26 

0.24 

0.27 

0.25 

0.28 

0.27 

Boizenburg 

48 

0.11 

0.17 

0.12 

0.19 

0.13 

0.20 

0.13 

0.21 

Braunschweig 

48 

0.19 

0.17 

0.18 

0.20 

0.26 

0.22 

0.20 

0.23 

Magdeburg 

48 

0.10 

0.17 

0.16 

0.19 

0.12 

0.20 

0.15 

0.21 

Wernigerode 

48 

0.08 

0.16 

0.10 

0.18 

0.10 

0.19 

0.10 

0.20 

Weissen 

48 

r  18 

0.19 

0.20 

0.21 

0.21 

0.22 

0.23 

0.23 

consistently  smaller 

for  the 

other 

three  stations 

The 

nesul ts 

;  for  46 

MVAR 

forecasts  for  the  third  transformed  variable,  visibility,  are  shown  in 


Table  7.  With  the  exception  of  those  found  for  Boizenburg,  the  RKSE's 
computed  from  the  forecasts  agree  ouite  well  with  those  expected  fro^  tree*-;/. 

Table  7.  Comparison  of  theoretical  forecast  RMSE's  (T)  with  those  computed 
from  actual  MVAR  forecasts  (C)  using  1976  transformed  visibility 
_ data. _ 

RMSE's 


No.  of 

3-hour 

FCST 

6-hour 

FCST 

9-hour 

FCST 

12-nour  :C 

Stati on 

rcrecasts 

C 

T 

C 

T 

r 

"r 

Hannover 

46 

0.09 

0.13 

0.09 

0. 15 

0. 18 

C.  16 

n  -  ~ 

V  .  -  - 

- 

Bremen 

46 

0.09 

0.14 

0.14 

0.16 

0.15 

0.15 

0.12 

r,  - 

Boizenb-rg 

45 

0.34 

0.16 

0.29 

0.18 

0.28 

0.19 

r\  7  t 

9  "7 

Braunschweig 

46 

0.10 

0.10 

0.09 

0.12 

0.21 

0.13 

0.  lr- 

3  " 

Maqdeburg 

46 

0.14 

0.11 

0.12 

0.14 

0.14 

0.16 

O  ^  9 

Wernioerode 

46 

0.25 

0.18 

0.19 

0.20 

0.19 

0 .20 

0.26 

^  • 

Wei ssen 

46 

0.12 

0. 13 

0.16 

0.15 

0.19 

0.16 

0.22 

0.  1 

It  was  found  that  during  the  forecast  periods,  a  much  larger  number  of  zero 
visibilities  (resulting  in  a  transformed  variable  value  of  one)  were  actually 
observed  at  Boizenburg  than  at  any  other  station.  In  Figure  2b  we  can  see 
that  sucn  observations  would  result  in  an  increase  in  variance  and  thus  in 
•-Up  pMSE’s  for  the  transformed  visibilities. 

The  RMSE  comparison  for  58  temperature  forecasts  is  displayed  in 
TaK1<"  p-  while  those  for  140  wind  forecasts  are  shown  in  Tables  9  and  10. 

We  car  see  in  Tables  8  and  9  that  the  computed  RMSE's  are  slightly  larger 
than  their  theoretical  counterparts  for  the  temperature  and  u-wind  MVAR 
forecasts.  In  no  case,  however,  are  these  differences  significant.  The 
v-wind  MVAR  forecasts,  whose  RMSE's  are  shown  in  Table  10,  are  consistently 


Table  8.  Comparison  of  theoretical  forecast  RMSE's  (7)  with  those  computed 
_ from  actual  MVAR  forecasts  (C)  using  1976  temperature  data. 

RMSE  1 


3-hour 

FCST 

6-hour 

fcst 

9-hour 

FCST 

12-hou 

r  CCS4 

Station 

Forecasts 

C 

T 

C 

T 

r 

T 

C 

f 

Hannover 

58 

1.65 

1.50 

2.25 

2 . 08 

2.94 

2.33 

2.79 

2.71 

Bremen 

58 

1.71 

1.53 

2.07 

2.02 

2.45 

2.31 

2.72 

2.57 

Boizenburg 

58 

1.61 

1.62 

1 . 98 

2.08 

2.31 

2.36 

2.25 

2.57 

Braunschweig 

58 

1.54 

1.42 

2.11 

1.93 

2.50 

2.24 

2.55 

2-53 

Magdeburg 

58 

1.74 

1.53 

2 . 29 

2. 00 

2.44 

2 . 

2.43 

2 .  i 

Wernigerode 

58 

2.18 

1.68 

2.53 

2.18 

2.82 

2.49 

2.81 

2.72 

Weissen 

58 

1.93 

1.79 

2.59 

2.28 

2.57 

2.56 

2.26 

O  “*  i ' 

better  than  those  expected  from  theory.  In  summary,  we  have  seen  in  Tables 
5-10  that  for  the  most  part,  when  tested  upon  1976  observations,  the  fore¬ 
casts  produced  by  the  six  ninth  order  variable-at-a-time  MVAR  models  (developed 
using  1972-1974  data)  agree  quite  well  with  what  one  would  expect  theoretically. 

Table  9.  Comparison  of  theoretical  forecast  RMSE's  (T)  with  those  computed 
_ from  actual  MVAR  forecasts  (C)  using  1976  u-wlnd  data. _ 

RMSE's  (msec"1) 


3-hour 

FCST 

6-hour 

FCST 

9-hour 

FCST 

12-hour  FCST 

Station 

Forecasts 

C 

r 

T 

C 

T 

C 

T 

Hannover 

140 

2.02 

1.58 

2.32 

2.09 

2.97 

2.42 

2.99 

2.71 

Bremen 

140 

2.56 

1.78 

2.87 

2.28 

3.26 

2.60 

3.42 

2.87 

Boizenburg 

140 

1.95 

1.96 

2.31 

2.23 

2.71 

2.49 

2.69 

2.66 

Braunschweig 

140 

1.74 

1.51 

2.31 

1.95 

2.74 

2.26 

3.01 

2.5C 

Magdeburg 

140 

1.52 

1.53 

2.34 

1.88 

2.65 

2.16 

2.87 

2.36 

Wernigerode 

140 

2.17 

1.91 

2.71 

2.27 

2.81 

2.53 

2 .  q  c, 

2.72 

Weissen 


140 


1.99  1.97  2.65  2.29  2.97  2.57 


2.77 


Table  10.  fonparison  of  theoretical  forecast  RMSE's  (T)  with  those  compute-: 
_ from  actual  MVAR  -forecasts  (C)  using  1976  v-wind  data. _ 

RMSE's  (msec-1) 


3-hour 

6-hour 

FCS~ 

9-hour 

rCFT 

2  -  ho 1 

u r  rC:8 

Stati on 

Forecast 

s  C 

T“ 

C 

T 

C 

T 

jr 

‘T 

Hannover 

i40 

1.32 

1.41 

1.41 

1.74 

1.78 

1.97 

1.80 

2.13 

Bremen 

1.33 

1 . 75 

2. 13 

2.17 

2.48 

2.41 

2.62 

2 . 62 

Bo  1  z ?  r'b*jro 

140 

1.29 

1.73 

1.45 

1.96 

1.59 

2.12 

1 . 68 

2.26 

Brauns;  Kwei'V 

140 

1 . 32 

1 .42 

1.43 

1.70 

1  7  7 

i  .  .*  / 

1.90 

1.72 

t.  .  .  .* 

Magdeourg 

>•  o 

1  or 

*.  •  L.  J 

1 . 30 

1.33 

*!  07 

1  •  J 

■y  r\  r' 

1  .  -•  J 

2.0i 

1  *,C 

t  •  .  ‘ 

Wernigerode 

140 

1.90 

1.93 

1.94 

2. 18 

2.25 

2.34 

7  nr. 

2.-7 

Meissen 

140 

1.45 

1.57 

1.56 

1.82 

1 . 38 

2.0: 

•  Oil 

L.  .  ^  '  j 

Fina11>,  in  the  next  table  we  will  demonstrate  how  confidence  interv 
can  ce  placed  about  MVAR  predictions  using  two  visibility  forecasts  -:a  ’e  ■■■''  • 
the  variab !e-at- a  -  time  mode1.  From  the  error  analysis  of  the  actual  forecas*  s . 
whose  re^vl uc  are  outlined  in  Tables  5  -  10,  we  conclude  that  in  practice  the 
M VAR  mode'  m&kt  predictions  mucn  ’ike  we  would  expect  from  theory.  Con^ide-cf 
intervals  (C.J.)  to  be  placed  about  the  MVAR  predictions  can  be  computed  by 
rvulpinly :;i?  a  constant  (varying  in  size  depending  on  the  percent  C.I.  desired) 
by  :-.ne  sgut'-'e  root  of  the  estimated  forecast-step  prediction  error  variance. 
r>’oee  th:  in  equivalent  to  multiplying  that  constant  by  the  theoretical 
?.M5w'  s  snowr  in  Tables  5-10,  we  can  see  that,  as  one  would  exDect,  larger 
arm  avger  C.I.’s  will  be  determined  as  the  forecast-step  is  increased. 

The  results  of  two  visibility  forecasts  made  using  the  variable- 
at-a-time  MVAR  model  are  shown  in  Table  11.  The  forecasts  and  confidence 
intervals  have  been  transformed  back  into  normal  units  (km).  The  12-hour 


26 


A 


Table  11.  A  comparison 
intervals  (in 
_  (A). 

of  MVAR  vi sibil ity 
parenthesis)  with 

forecasts  (F)  and  800  confidence 
the  actual  observed  visioilities 

Forecast  Time  -  187. ,  March  29, 

1976 

3-hour 

6-hour 

9- hour 

12-hour 

Station 

F  A 

F  A 

JL  _A_ 

J_  _A_ 

20+  20+ 

7.8  20 

7.1  8 

5.3  7 

Hannover 

(3.9,20+) 

(3.1,20+) 

(3.0,20+) 

(2.5,20+) 

20+  20 

20+  20 

11  20 

5.7  8 

Bremen 

(3.7,20+) 

(3.2,20+) 

(3.1,20+) 

(2.6,20+) 

7.5  18 

8.2  15 

6.7  12 

4.5  6 

Boi zenburg 

(3.0,20+) 

(2.8,20+) 

(2.6,20+) 

:  5.1,20+) 

20+  20 

20+  15 

9  15 

r>  r\ 

0.0  Z? 

Braunschweig 

(4.3,20+) 

(3.3,20+) 

(3.5,20+) 

(3.3,20+) 

20+  16 

20+  16 

6.9  12 

6.5  12 

Magdeburg 

(4.1,20+) 

(3.5,20+) 

(2.9,20+) 

'  ;•  7  ?oi } 

V  4  •  •  **-  «.' 

4  20 

0.5  0 

0.6  0 

4.9  12 

Wernigerode 

(2.0,20+) 

(0.,  1.3) 

(0.,  1.5) 

(2. 1 ,20 _ ) 

20+  20 

7.3  20 

5.5  15 

5.7  6 

Weissen 

(3.6,20+) 

(3.0,20+) 

(2.6,20+) 

(2.6,20+) 

Forecast  Time  -  i 

06Z ,  April  16, 

1976 

3-hour 

6-hour 

9-hour 

12-hour 

Station 

F  A 

F  A 

F  A 

F  A 

1.1  0.4 

2.1  6 

3  10 

3.2  8 

Hannover 

(0.6, 1.7) 

(1.3, 3. 8) 

(1. 7,7.7) 

(1.8,20+) 

1.3  1.5 

2.2  6 

2.5  8 

3.1  9 

Bremen 

(0.7,2. 1) 

(1.3, 4. 2) 

(1.4, 5. 2) 

(1.7,20+) 

0.7  0 

1.4  7 

1.8  0 

2.3  10 

Boizenburg 

(0.2, 1.4) 

(0.6, 2. 6) 

(0.9, 3. 7) 

(1.1, 5. 5) 

0.9  0.1 

1.9  6 

2.3  7 

3.0  6 

Braunschweig 

(0.5,1. 3) 

(1.2, 2. 9) 

(1.5, 3. 8) 

(1.8, 6. 2) 

0.7  3.5 

1.5  6 

1.8  8 

2.1  6 

Magdeburg 

(0.3, 1.2) 

(0.8, 2. 4) 

(1.0, 3. 2) 

(1.2,4. 1) 

1.8  4.5 

1.8  4 

3.1  6 

3.2  6 

Wernigerode 

(0.9, 3. 4) 

(0.8, 3. 8) 

(1.5,20+) 

(1.5,20+) 

1.0  4 

1.8  6 

2.2  10 

2.6  8 

Weissen 

(0.5, 1.7) 

(1.0, 3. 2) 

(1.3,4. 3) 

(1.4,6. 1) 

period  after  18Z  on  March  29,  1976,  was  basically  a  time  o  +  high  visik’'' 
In  this  case  ev>-ry  80'  C.I.  about  the  MVAR  forecasts  included  the  actual 
visibility.  It  is  especially  notable  that  the  model  predicted  very  we1’ 
the  very  low  visibilities  actually  observed  at  Wernigerode  6  and  9  hourr 
a'ter  tne  initial  time.  Tne  12-hour  period  following  06Z  on  April  16, 
i  /i ,  was  one  in  which  visibilities  were  very  '  ow  after  3  hours  •  •  t.r 
steadily  improveu  over  the  rest  u+  the  period.  The  very  low  visit'’  ‘t'i-  • 
at  the  3-‘>0!:r  r-r  k  were  we"1!  forecast  and  for  every  station  fnp  vis't  i  . 
were  predicted  to  improve  out  tc  12  hours.  In  th’s  case  the  model  c;:  n<v 
improve  the  visibi1 ities  as  fast  as  nature  and  only  a  few  or  the  3C  C.I. 
include  the  actual  observations.  prom  this  table  we  can  see  tna*  the 
confidence  intervals  can  probably  be  best  used  by  an  actus’  forecasts-  t: 
specify  a  minimum  expected  visibi ’ity  in  high  visibility  situations  and  • 
m..:< '■  i-m  expected  visibility  in  "low  visibility  situations.  This  type  o* 
interpretation  of  the  C.I.'s  is  a’ so  appropriate  for  the  oth?"  two  trans¬ 
formed  yciri^'.-’es  (ceiling  height  .  nd  height  of  first  cloud  layer),  "he 
customary  interpretation  of  the  C.I.'s  as  a  range  in  which  we  expect  the 

variable  to  lie  can  be  applied  to  the  other  three  variables  (temperature, 

,  o 

u-wind,  and  v-wir.d).  This  range  was  found  to  be  the  order  of  +2  C  and 

+2  m/s:c  ’or  a  o-hour  forecast  of  temperature  and  the  wind  components, 

'ecpect i v  ’y.  For  12-hour  forecasts  it  was  found  that  the  temperature  and 

v:  r  ’  ranges  were  about  +  3°C  and  +3.5  m/sec,  respectively. 
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5.  Conclusions 


In  this  study  multivariate  autoregressive  climatological  models 
were  developed  to  be  used  for  short-term  forecasting  of  six  atmospheric 
variables  (temperature,  visibility,  u-wind,  v-winb,  ceiling  height,  and 
height  of  the  first  cloud  layer)  over  a  limited  area  in  a  tactical  situation. 
After  a  preliminary  data  analysis  it  was  found  that  the  cloud-height  variables 
and  visibility  could  best  be  forecast  if  they  were  first  transformed  so  that 
the  variance  of  low  ceiling  and  visibility  situations  be  emphasized  over 
that  of  high  ceiling  and  visibility  conditions.  Various  forecast  models 
were  investigated,  and  it  was  found  that,  given  the  tactical  constraints, 
the  best  models  were  ninth  order  variable-at-a-time  MVAR  models  in  which 
an  observation  vector  consisted  of  the  values  of  the  variable  to  be  predicted 
at  the  seven  German  stations  used  in  this  study  (Hannover,  Bremen,  Braun¬ 
schweig,  Boizenburg,  Magdeburg,  Wernigerode,  and  Weissen).  The  oredicticn 
error  variance  matrix  and  the  coefficient  matrices  for  each  of  the  six 
variable-at-a-time  MVAR  forecast  models  are  tabulated  in  Appendix  A. 

Using  these  models  one  can  expect  to  make  3-hour  forecasts  which  explain 
approximately  95%,  75%,  60%,  55%,  45%  and  45%  of  the  variance,  respectively, 
for  temperature,  u-wind,  v-wind,  visibility,  ceiling  height,  and  height 
of  the  first  cloud  layer. 

The  MVAR  models  were  developed  using  data  collected  during  1972- 
1974.  The  models  were  then  tested  independently  using  1976  observations 
and  it  was  found  that  the  actual  forecast  errors  agree  quite  well  with 
what  would  be  theoretically  predicted.  Using  the  estimated  prediction 
error  variances,  confidence  intervals  were  determined  to  be  placed  about 
MVAR  forecasts.  It  was  found  that  80%  C.I.'s  of  +  2°C  and  +  2  m/sec  cou’d 


be  placed  about  the  3-hour  forecasts  of  temperature  and  the  wind  component 
respectively.  Due  to  the  variable  transformation  made  upon  the  visibility 
and  cloud  height  variables,  it  was  found  that  the  confidence  intervals  cou 
be  best  used  to  determine  minimum  expected  visibilities  (or  cloud  heights 
in  high  visibility  (or  ceiling)  situations  and  to  determine  maxi mums  in  th 
poor  visibility  (low  ceiling)  situations. 
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COVARIANCE  AND  COEFFICIENT  MATRICES  FOR  THE  SIX 
VARIABLE-A7-A-TIMC  MVAR  MODELS 

VARIA8LE-AT-A-TIHE  MVAR  MODEL  FOR  T’ttwoeRATUKE  MODEL  ORDER  IS 


ONE-STEP  PREDICTION  ERROR  COVARIANCE  MATRIX 


2. 

2390 

1. 3424 

1  .  1723 

1. 3838 

1.0775 

1.4253 

1.2331 

1 . 

3424 
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0.5150 

-0. 2799 
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0 

.  e  Fan 

-0.0041 

-0.0218 

-0.  0147 

C.  0  039 

-C. 0099 

0.0  176 
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COEFFICIENT  MATRIX  A5 
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■'A3IABLE-AT-A-TII1E  MVAR  MODEL  FOR  U-rfINC  MODEL  ORDER  IS  9 


ONE-STEP  PREDICTION  ERROR  COVARIANCE  MATFIX 
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j  • 

r  +  p ;  C 

0.1467 

-C.  27  2  0 

-0.2264 

-0.0510 

-0.0798 

-0.0474 

-o. 

9  ►- 

0 .2544 

-0.2311 

-0.0330 

-0.3417 

-0.0463 

-C. 1305 

_  r\ 

0  14  1 

0.2349 

-0.  1071 

-0.0387 

-0. 1023 

-0. 3C61 

-0. 1609 

-0. 

0  5  4  3 

0.2357 

-0.  1238 

-0.0290 

-0. 0479 

-0.0932 

-  0.  4301 

-0. 

5  129 

0. 1582 

-0. 2172 

-0.1432 

-0. 0874 

-0.1719 

1 

o 

o 

& 

o 

-0 . 

2  4  5-7 

COEFFICIENT  MATE 

IX  A 2 

0.008C 

0.0581 

0.0 180 

0. 0174 

C. 0125 

0.  05Q5 

0.  0  19  9 

9 . 0  4  9  1 

-7. 0430 

0.02C5 

0.0535 

0.0471 

0.0314 

0.  0047 

0.0075 

0.0120 

-0.  1045 

-  0.  0023 

0.0619 

C.  0 177 

0. 0094 

0.061? 

0  .  106  4 

-0 . 0077 

-0.  1098 

0.0066 

C. 0564 

-0.0049 

0.0414 

0.0627 

-0.0246 

0.0138 

-0.0476 

C. 0152 

-0.  0005 

0  0«6  6 

o. 0273 

-O.01H 

0 . 0006 

0.0458 

-  0.  083  9 

0. 021  9 

0.O821 

0 . 0296 

-0.0358 

-  0.  Cl  26 

-0.C  1 16 

0.  069  1 

-0. 0518 

f-FFTCIENT  MATRIX  A3 


C .: U  7  7 

-0 . 0020 

0.0091 

0. 0094 

-0.0027 

-0.0165 

-0.  CO  13 

C  .  0  2  9  1 

-0.0370 

0.0251 

0.  0404 

0.0  112 

-0.  0030 

-0.0153 

0  .  C  2  6  9 

0.0130 

-0.0481 

0. 0176 

-0.0011 

0.0569 

-0.  0  250 

0 . 0754 

-0.0046 

0.0231 

-0.  0246 

0.0279 

-0. 0066 

0.C  1  1  8 

-O.c 1 32 

0.0307 

0.0182 

0. 0715 

-0.0062 

-0.0065 

-0.  0021 

0.066  ’ 

-0.0107 

-0.0140 

0.  0422 

0.  0350 

-  0.  0608 

0.0077 

-  0  .  C  2  7  i 

0.0237 

0.0172 

0.0324 

0.0597 

0.0106 

-0.031  1 

-  r FT ICIENT  MATRIX  AU 


-C .0268 

-0.0154 

0.0293 

-0.  0299 

0.0155 

0.  0369 

0.3237 

-0.0157 

-0.0432 

0.0154 

0.0103 

-0.0023 

0. 003  1 

0.0075 

0.0091 

-0,0106 

-0.0490 

0. 0135 

-0.0114 

-C. 0088 

0.0 1 35 

-0.0254 

C . 00  06 

0.0221 

-0.  0672 

0. 0129 

0.0175 

-O.C06C 

-0.0187 

0.0006 

-0.0145 

-0.  0039 

-0.0218 

-0. 0020 

-0.  3032 

-0.0165 

-0.0304 

0.0464 

-0. 0020 

-0.0038 

-  0.  0007 

0.C075 

0.0082 

0.0172 

0.0341 

-0.0046 

-0.1643 

-0.0032 

-0.  0  369 

34 


COEFFICIENT  MATRIX  A5 


-0 . C44fl 

0.0321 

-C. 0225 

0. 0272 

-0. 0C85 

-0. 0072 

-o.om 

-0.0182 

0.0039 

0.0045 

0.  Cl  99 

0.0C  13 

0. 0004 

-0.CO39 

-0.C468 

0.0395 

0. C 144 

-0. 0129 

0.0461 

0.0196 

0.0215 

0.0215 

0.0375 

-0.0011 

-0.0141 

0.C321 

-0.0114 

0.  0  225 

-0.0147 

0.0142 

0.0173 

C. 0547 

-0. 0  115 

-0.0199 

0.  Cl  88 

0.0223 

-0.0162 

0.0144 

0. 02  19 

-0.0109 

-0.0363 

0.  0C72 

0.0039 

0.0202 

0.C155 

0.0376 

0.0638 

-0.0046 

-0.0167 

COEFFICIENT  MATE 

IX  A  6 

0.0055 

-0.  0258 

0.0190 

-0. 0018 

-C. 0  172 

-C. C085 

-  0.  0046 

-0.0140 

-0.0432 

0.0188 

-0. 0231 

-0.0050 

0.0017 

0.  0080 

0.0589 

0.0006 

-0.028  1 

-0. 0538 

-0. 0589 

-0.0130 

0. 0125 

0.0228 

-0.0315 

0.0269 

-0.  0719 

-0.0252 

-C. 0009 

-C.  00  34 

0.0381 

0.0182 

-0.0098 

-0. 0497 

-0.0267 

0.0174 

-0.  0141 

0.0399 

0.0072 

-0.0049 

-0. 01  13 

-0.0572 

-C . 0224 

0. 0030 

0.0228 

-0.01 1C 

-0.0C1  1 

-0.0349 

-C.0219 

0.0214 

-0. 0224 

COEFFICIENT  MAT: 

FIX  A  7 

-0  .C298 

C.  0086 

-C.0111 

-0.0251 

0.0295 

-  0. C 080 

-0.  C244 

-0 .0 162 

-0.C088 

-0.0249 

-0.0068 

C .  0548 

-0.017  1 

-0.0347 

—  ^  •  G  4  8  9 

0.0186 

-0.0437 

0. 0261 

0.0083 

C. 0150 

-0.  0074 

-0.0105 

0.0144 

-0.0115 

-0.  0448 

0.0281 

0. 0318 

-0. 0120 

-0.0763 

0.0255 

-0.0211 

0. 0390 

C . 0057 

-0.0071 

0.  0006 

0.0007 

-0.0319 

-0.0396 

0. 0076 

0.0229 

-0.008  1 

-0.  0  1  19 

-0.0496 

0.0135 

-0.0132 

-0. 0086 

0.0305 

-0.0133 

-0.C208 

COEFFICIENT  MATRIX  A8 

0.0021 

0.  0168 

-0.0C94 

-0. 0379 

-0.0445 

-0.  0104 

-0.  C001 

0.0097 

-0.0021 

-0.0125 

-0. 0654 

-0.0207 

0. 0045 

-o.oooc 

0.0075 

0.0203 

-0.0400 

-0. 0348 

-0.0200 

-0. 0C5C 

-0. 01 34 

-0.0246 

-0.0022 

-0.0419 

-0.0412 

-0.0549 

-0.0281 

-0.  0163 

0.0001 

0.0049 

-0.0082 

-0.0233 

-C. 0463 

-0.0154 

-0.C171 

0.0020 

0.01  19 

-0.0180 

-0.  0249 

-0.0120 

-0. 0521 

0.0042 

-0.0064 

0.0272 

-0.0120 

-0. 0476 

-0.0252 

0. 0106 

-0. 0236 

COEFFICIENT  MATRIX 

A9 

-0.C081 

0.0059  - 

0.01C3 

0. 0889 

0. 0334 

0. 0131 

-0. C060 

C.C318 

-0.0134  - 

0.0012 

0. 0670 

0.0375 

-0.  01C5 

-0.  00  36 

0.0088 

0.0152 

0.0060 

0. 0729 

0.0  154 

0. 0097 

0.  C076 

0.0161 

0 . 0444 

0.0052 

0. 0427 

0.0366 

0.0  375 

0. 0094 

-0.0015 

0.0165  - 

0.0093 

0. 0781 

0.0368 

0.0154 

0.  0058 

-0.0197 

0.0174  - 

0 . 0024 

0.0786 

0.0335 

-0.0125 

0. C 194 

0.0198 

0.0C32 

O.C  103 

0. 0458 

0.0264 

C. 0465 

-0.0012 

VARJA3LE-AT- A-TIME  MVAP  MODEL  FOR  V-WINE  MODEL  ORDER  IS  7 


ONE-STEP  PREDICTION  ERROR  COVARIANCE  MATRIX 


2 . 0  0  Oil 

0.395C 

0.5895 

0. 7087 

0.5176 

0.5277 

0 . 

45  2  r 

C . 89  55 

3.0762 

0.5891 

0.4587 

0.4061 

0.3197 

0. 

4  3  3  7 

0.5395 

0.589  1 

2.9894 

0. 5059 

0.4896 

0.2117 

0. 

5  84  2 

C .7087 

0.4587 

0.5059 

2. 0322 

0. 6705 

0.5179 

9. 

567  2 

0.5176 

0.4081 

0.4896 

C. 6705 

2.5  1  16 

0.6577 

0. 

f  r\n  c 

0.5277 

0.3197 

0.2117 

C  .  5 1 73 

0.6577 

3.7166 

n  p 

1*  ^  R  1 

0 . 4520 

0.4338 

C  .5342 

0.56 7 J 

0.6CC9 

C  .  4  3  5  8 

u  •'  "  •' 

CO  EE FTC 

I  ENT  MAT? 

TX  A1 

-0.3342 

-0.  3204 

-0 . 0893 

-0.  0854 

-0.0345 

-0.0  29  1 

^  »■  ,  ..  . 

-0 . 1 306 

-0.65  86 

-0.0335 

-C  .  04  9o 

-0.0284 

0.016? 

-0.  •  i  4  ■ 

-0.  1214 

-0.2677 

-0.2575 

-0.0820 

-0.0660 

-0 . 022° 

—  P,  r  •**  ^  ^ 

-0.2151 

-0.139  1 

-0 . CQ4 1 

-0. 2721 

-0.0707 

-0.0367 

-  \  0  8 

-0.  1 6  4  C 

-0.1425 

-  0.  0833 

-0.1899 

-0.2004 

-  C  .  0  5  3  7 

—  ^  «  -*  -• 

-0.1973 

-0. 1506 

-0.0778 

-C. 0992 

-0.0671 

-  C.  3 54 >> 

-  #  ^  :  '0 

-0. 1621 

-0. 1806 

-0.1 309 

- : . 1272 

-0.0732 

-C.04E5 

- " .  ;v  r 

COEFFICIENT  MATE 

IX  A  2 

-C.0056 

0.0325 

C.C315 

-  C  .  03  4  0 

0.0065 

-0. CC07 

0.  005  2 

C .0366 

-0.0160 

C.C27B 

C. 0284 

-0.0067 

-0.0059 

0.r37i 

0.0833 

0.0584 

-0.0877 

-0.  03  42 

-0.0100 

0.0007 

0.  0  24  4 

0.026  1 

0.0579 

O.C3C8 

-0. 0791 

0.0007 

0. 0009 

n .  0  r.  1  (:) 

C.0104 

0.  0682 

0.0152 

0.0018 

-0.0284 

0.0214 

0.  0  0»7 

-C. 0228 

0. 0656 

C .0278 

".0783 

0.0275 

-  0.  0638 

-  0 . 0  ■  2  ’ 

0  .  C  4  0  9 

0.0357 

0.0153 

-  .0127 

-0.0096 

C  .  0  2  1 6 

-  0  .  0  4  2  8 

COEFFICI 

ENT  MATRIX 

A3 

-0.0017 

0.0574 

0.0284 

-0.  0081 

0.0183 

-0.0007 

-0.0162 

-0.0064 

-C.0177 

0.0256 

0.0170 

0. 0293 

0. 0099 

-r,.  Cl  28 

0.00 9 9 

0.0302  - 

0.0459 

0.0005 

C . 0336 

0.0287 

-0.  ODBC 

0.0301 

0.  0472 

0.0306 

-0. 02  12 

0.0224 

0.0127 

0.  004  1 

0.0225 

0.0198 

0.0162 

-0. CO  85 

-0.0404 

0.0280 

-0.0042 

-0.0034 

0.0453 

0.0034 

0.0218 

-0.0 167 

-0.0177 

0. 0 1 84 

0.0242 

0.0256  - 

0.0039 

-0. 0132 

0. 0096 

0.0174 

-0.  0  08  5 

COEFFICIENT  MATR 

IX  A 4 

-0.0236 

0.0071 

0.0011 

-0. 0063 

-0.0 128 

-C. 0152 

0.0044 

-0.0106 

-0.0060 

0.0270 

0.0038 

-0. 0036 

0. 0033 

0. 01 38 

0.0025 

0.0100 

-0.0420 

-0.01  09 

-0.0117 

0.0134 

0.  0 155 

-0.0129 

0.0040 

-0.0203 

-0.0313 

-0. 0170 

-0.0091 

-  C.  0  1  09 

0.0101 

-0.0116 

0.0174 

-0.0158 

-0.0042 

0.0010 

-0.  000  3 

-0.0053 

-0.0448 

0.0275 

0. 0373 

-0.0094 

-0.0129 

0. 0016 

0.0114 

0. 0020 

-0.0011 

-0. 0099 

-0.0115 

-0.  0026 

-0.  0056 

36 

COEFFICIENT  MATRIX  A5 

1 

-3.0069 

0.0035 

0.0215 

-0. 01 36 

0.0115 

C. 0020 

0.0239  1 

0.0104 

-0.0336 

0.0057 

-0. 0020 

0.0125 

0. 0036 

9.0101  |J 

-0.0121 

0.0314 

-0.0472 

0. 00  16 

0.0090 

0.0331 

0.  0097  fl 

0.0043 

0.0057 

0.0260 

0.0030 

0.0327 

C. 0286 

-0.  0238  .■ 

-0.0127 

-0.0180 

0. 0186 

-0. 0089 

-0.0064 

0.0130 

-C.O 161  {■ 

-G .0280 

-0.0064 

0.0028 

-0.0157 

0.0343 

0.0063 

-0.0290 

Q . 00  00  -0.0069 

COEFFICIENT  MATR 

0.0166 

IX  A6 

0.0123 

0.0023 

0.0243 

-0.0068 

-G.C037 

-0.0127 

0.0210 

-0. 0232 

C.O  106 

-0. 0079 

0.0159  1 

-0.0167 

-0. 0399 

-0.0036 

-0.0265 

0.0224 

0. 01 30 

0.0067  I 

C .00  16 

0.0063 

-0.0415 

-0. 0043 

C. 0 172 

-0. 0049 

0.0214  I 

0.0C52 

-0.0124 

-0.0179 

-0. 0345 

0.0019 

C. 0057 

-0.  000  2  Cl 

0.0076 

0.0060 

-C.01 47 

-0. 0128 

0.0013 

-0. CC7  0 

0.  0353  I 

-0.0137 

0.012C 

0.0088 

-0.0069 

0.0  189 

-0.0132 

0.049C  1 

0.0098  -0o  016  9  -0.0187 

COEFFICIENT  MATRIX  A7 

-0.0255 

0.0085 

-  0.  0005 

0.0129  I 

-C  .0290 

-0.0063 

-0.0177 

-0. 0269 

0.0C56 

0.0090 

-0.0148  1 

-0.0026 

-0.0264 

-0.0151 

-0.  03  11 

C. 0  1  14 

0.0030 

-0.  0080  1 

-C  .0196 

-0.0179 

-0.0245 

-0.0351 

-0.0C33 

-0.0187 

-0.0085  1 

-0.0245 

0. 0056 

-0.0144 

-0.02  97 

0.0023 

-0. 0079 

0.  0002  ] 

-0.0181 

0.0070 

-0.0062 

-0.0181 

-0.0282 

-0.0104 

-0.  0052  1 

-0.0016 

-0.0036 

-0.0085 

-0. 06  06 

0.0279 

-0. 0553 

-0.0498  1 

-0.0175  -0.0083  -C.0075 

COEFFICIENT  MATRIX  A8 

0. 0028 

-o.oo eo 

-0.0090 

-0.01 17 

1  , 

-C .0033 

0. 0189 

-0.C198 

-0. 0384 

-0. 0166 

-0,0085 

0.0097 

-0.0050 

-0.0459 

0.C1 17 

-0.  0132 

-0.0241 

0.0067 

0.  0039 

0.0428 

0.0085 

-0.0308 

-0.0353 

-0.0390 

0.0027 

-  0.  0  067 

-0.0046 

0.0146 

-0.0358 

-C.  04  09 

-0. 0354 

-0. C  164 

-0.0018 

-0.0050 

0.0167 

-0.0262 

-0. 00  15 

-0.0292 

0.0007 

-0.0032 

C.0681 

0. 0169 

-0.0238 

-0.  0095 

0.0130 

-0.0763 

0.0063 

0.0155  0.0087  -0.0180 

COEFFICIENT  MATRIX  A9 

-0.  0325 

-0.0276 

-C. 0139 

-0.0061  j 

ij 

0.0251 

0.0302 

0.0096 

-0.0255 

O.C  128 

-0.0001 

0.0178 

0.0359 

0.0198 

0.0058 

-0.0274 

0.0127 

-0.0091 

0.0255 

0.0644 

0.0422 

-0.0255 

0.  0033 

0.0041 

-0.  0008 

0.0191  j 

0.0357 

0.0472 

0.0195 

0.0156 

0.0039 

-0.0078 

O.C 146  I 

0.0250 

0.0454 

-0.C018 

-0. 0074 

C.O 1C8 

-0. 0055 

0.0110  Jft 

C.0162 

0.0448 

0.0169 

0.0025 

-0. C  187 

-0.0113 

0.  0225  M 

0.0367 

. 

0.0289 

0.0098 

0.0131  0.0121 

37 

0. 0034 

1 

o  o  o 


■•.nT^Tr.jT-j.Tpr.  1VAP  5DEL  FOR  VISIBILITY 


MODEL  OPT*'? 


?*’  F-STEP  PREDICT  TON  ERROR  COVARIANCE  MATRIX 


*  r  i  j:  -? 

w*  <  •  <  V  < 

C.GC5C 

9.CC33 

0.0052 

0.0027 

C  .  C  0  1  6 

0.  CC30 

0.0C5C 

0.0197 

0.0035 

0. 0023 

0.0016 

0.0011 

0. ^024 

0 .0033 

0.0035 

0.0241 

0. 0024 

C.004C 

0.0012 

0.0065 

0.0052 

0.0C23 

C. 0024 

0. 0097 

0.0030 

0. 00  18 

0.00  14 

•\  r  •>  *■» 

vy  *  •'  J  <L  . 

C .0016 

0 .0040 

0.0030 

0. 0 132 

C. 0013 

0.  r  0  1 . 

0.0016 

0.0011 

0. 0012 

0.0013 

0.0018 

C  .  0  3  3  5 

0 .  r  o  r  ' 

'•  A  O  T  ^ 

*  ’ 

.  00  ->4 

9.0C65 

3.0019 

0.C034 

0 .  0  0  :  ? 

O  /\  «  •*» 

CCEPFIC.7. 

EF'T  ’■*  a "  _r - 

V  }  ’ 

-0.4523 

-0  .  179  5 

-0.0362 

-0. 1164 

-0.G43C 

-0.033° 

_  'A  ' 

-  0 .  17  9  5 

-  ■> .  u  - 23 

-  0.02  9  3 

C. 01 99 

-0.0325 

-  c .  c  i  e  3 

_  C-  r  ~  '  • 

-0.0864 

■0.  139  6 

-C . 4253 

-0.0256 

-0.0959 

-  0.  005  1 

_  r i  *  t  n  r. 

-C. 71 51 

-0.1018 

-  0. 0162 

-0. 3845 

-0. 0901 

-0.0316 

-o.  o  i6f- 

-0. 1393 

-0 .0449 

-0 . 0438 

-0. 0984 

-0. 5623 

-  0  .  0244 

_  '3  -  r  .■ 

-0.0476 

0.0182 

0. 01  Ci 

-0.  1014 

-0. 1 071 

-0. 4005 

0.  C  18  0 

-0 . 0497 

-0.090C 

-0.1273 

-C . 0033 

-C.  0  8  54 

0. 0034 

s'-  ' 

jprpt1  Vi 

01  EFT  MATRIX 

A2 

0.03  1? 

0 .0495 

3.0288 

0.0392 

-0. 0098 

-0.0028 

7.9  32  ’ 

.  o 1  o  6 

0.02  3  5 

0.013 8 

0197 

-0.0234 

-C . 0C4  1 

C  22  1 

0.10c’ 

,•>  A’.  *  1  _ 

0.01 07 

0.  04  22 

C. 0034 

-  0. 0 0° 2 

0.  0  3°2 

0.0  32  1 

0.0600 

0.0097 

9. 02  63 

0.0087 

-0.0043 

0 .  C  ’  8  2 

0.058  1 

9.9546 

9.0238 

2.0501 

-C .  0  627 

-0.0147 

0.  C  2  ? i 

3.  0  2  .32 
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V  • 

0  0  2  5 

-0. 

0  269 

0.  0  366 

.  r.  ■;  r; 

0 . 0  '  T  -  c . 

0  0  9  1 

X/  •  •  'vJ  ^  () 

-0. 

C  i  92 

0. 

0  279 

-  0 .  C  3  6  3 

•  97.- ?'  C 

I  7  NT  K  AC  P 

IX  A  9 

0.0394 

0 .0110 

-C.0C6? 

0. 0022 

-0.0  120 

0.0018 

0.  0090 

0.0941 

0.  00  4  c 

0.001  1 

-0.0211 

-  3  .  0  2  6  9 

0.0281 

0.0  138 

'  .0138 

-0 . 00  48 

-0.080R 

0.  01.32 

-0.0188 

0.0160 

-0.0182 

'.0  216 

0. 00  14 

-0.0014 

3.0146 

-0.0342 

-0. 02  18 

0.  0093 

0.0'.  30 

•0.00  31 

-C  .0130 

-0.  C J 32 

-C. 0 433 

-0.0324 

0.  0.037 

0 . 0169 

-0,0347 

0.0293 

-0. 0167 

0.0030 

-0.0012 

-0.  0  298 

0.0 .32 

0.00  78 

-0.0807 

C. 0498 

0.0  140 

-  0.  0383 

-0.0080 

COEFFICIENT  MATRIX  A5 


-0.0102 

-0.0198 

0.0541 

0. 0081 

-0.0216 

0. 0146 

-0. 0086 

-0.0012 

-0.0223 

0.0189 

0. 0209 

-0. 0395 

-0. 0093 

-0.  0068 

0.0302 

0.0012 

0.0351 

-0.  0275 

-0.0175 

0. 0386 

0.  00  7 u 

-0.0093 

-0.00  20 

0.0459 

0. 0180 

-0.0090 

0. 0  0°  4 

-0.0111 

-0.0200 

-0.0002 

0.0282 

0. 02  63 

-0.0171 

0.0131 

0.  00  35 

0.0223 

0.0112 

0. C356 

-0. 0451 

-0.0540 

0.0019 

0.  0564 

-0.0058 

-0.0026 

0.0608 

-0. 0198 

-0.0201 

0. 0009 

0.  0  070 

COEFFICIENT  MATRIX 

A6 

-0.0078 

-0.0202 

0.0131 

-;).  00  88 

0.01C6 

0.0134 

-0. C  255 

0.0422 

-0.0241  - 

0.0163 

-0. 0453 

C. 0242 

0. 0336 

0.0126 

-0.0139 

-0.0126  - 

0.0146 

0. 01  15 

0.0010 

-0.0112 

-0. CC44 

-C.C041 

-0. 0253 

0.0230 

-0. 0046 

C . 0116 

-0.  0016 

-0.  01 12 

-0.0090 

-0.  009  3 

0.0140 

0.  03  81 

-0.0171 

-0.0098 

-C. C 230 

-0.01 08 

0.0C39  - 

C.0089 

0. 0209 

0.0221 

0. 0183 

-0. 0330 

0.0068 

-0.0110 

0.0152 

0.07  31 

0.0051 

0.0268 

-0.  0  197 

COEFFICI 

ENT  MATRIX 

A  7 

-0.0102 

0.0080 

0.0017 

0.0015 

0.0079 

0.0046 

-0.C257 

0.0061 

-0. 0070 

0.0233 

-C. 0032 

-0. C  163 

0.0026 

-  0 .  C  4  4  0 

-0.0159 

0.0047  - 

0.0115 

0. 0073 

-C.CC09 

-C.C155 

-0.0099 

0.0001 

-0.0076 

0.0234 

-0.  0291 

-0.0151 

0.0392 

-0. 0025 

-0.0159 

-0.0001 

0.0258 

-0.  0070 

0.0092 

0.0205 

0.0219 

-0.0099 

-0. 0085 

0.0005 

-0. 0094 

0.0016 

-0.0398 

0.0144 

-0.0179 

-0.0041 

0.0082 

-0.0250 

-0.0023 

0. 0071 

0.0029 

COEFFICIENT  MATRIX  A8 

-0.0385 

-0.04  06 

0.0325 

-0. 0443 

-0.0326 

-0.0113 

-0.  C  187 

-0.0352 

-0.0563 

0.0052 

-0. 0354 

-C.0C22 

-0. 0284 

0.0094 

0.0233 

0.0032 

-0.0255 

-0.  0208 

-0. 0088 

-0.0351 

-0.  0090 

-0.0230 

-0.0087 

0.0009 

-0.0382 

-0.0332 

0. 0036 

-0.0443 

0.0123 

-0.0088 

-0.0022 

-0. 0015 

-0.0  7  34 

-0.0240 

-0. 0092 

-0.0240 

0.0090 

-0.0066 

0.0121 

0.C012 

-0.  0307 

C.C107 

0.00  15 

-0.0104 

0.0316 

-0. 0524 

-0. G 154 

0.0251 

-0.  0407 

COEFFICIENT  MATRIX  A9 

0.0324 

0.0342 

-0.0214 

0.  0381 

-0.0010 

-C. 0090 

0. 0240 

0.0149 

0.0318 

0.0003 

0. 0285 

-0. 0204 

0. 0126 

0.0125 

-0.0125 

0. 0038 

0.0088 

0. 0194 

0.0046 

0.0366 

0. 0027 

0.0226 

0.0088 

0.0297 

0.  0280 

0.0  1 18 

-0. 0046 

0.0151 

-0.0144 

0.0276 

-0.0103 

0. 0033 

-0.0231 

0.0097 

-0.  0032 

0.0011 

0.0291 

-0.0030 

0.0183 

0.0166 

-0.0314 

-0.  0127 

-0.0106 

0.0162 

-0.0089 

0.0411 

0.0356 

-0.0249 

0.  006  1 

or  -;nix  '■ 
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T”:s  P  a  0  j  H  A  M  IS  USED  TO  SELECT  ’’HE  PROPER  ORDER  HVAE  MODEL  AND 
LET  ERRIN  ES  THE  COVARIANCE  AND  COEFFICIENT  MATRICES  FOR  THAI  MODEL 
GIVEN  II  NP-DI MSNSICNAL  OBSERVATION  VECTORS.  THE  MAXIMUM  ORDER 
MODEL  CONSIDERED  IS  LG  OLD- 1 .  AFTER  THE  FINAL  MVAR  MODEL  HAS 
SEEN  DETERMINED,  THE  PREDICTION  ERROR  MATRICES  FOR  V ARI 0  US -ST E? 
FORECASTS  ARE  COMPUTED  BY  SUBROUTINE  ERRVAR. 

n:  MEN3IH-J  GAM  (7, ->,10)  ,A  (7,7. 10)  .A3  (7,7,  10)  ,B  (7,7,  1C!  , 

I  B”  (  7,  '•  .  . .  T  V  7,  976  6)  ,  X  M  N  '">)  ,  S  (7,7)  ,33  (7,7)  ,  EE  (7,7)  ,  .T:  (7 , 7 

/  0(7,7'  ;  1  '(*  0)  ,  .*■  A  (?  *7}  ,  HO"  (7, 7) 

np -  ;> 

Np  IS  V*  M'-Rf?  0?  PREDICTOR  ’  R  I  A  3L  ES . 


*<<.  : 


*IC  IS  THE  NUMBER  0,?  VAR  I  A?L~?  TO  BE  PREDICTED,  NC  II  LESS 
EQUAL  TO  NP. 

LGOL  D=  1  0 

LG  OLD  IS  ONE  MORE  THAN  THE  MAXIMUM  ORDER  MVAP  MODEL  r'D  EE 
LG=LGOL  D 
11=8768 


II  IS  THE  NUMBER  OF  OBSERVATION  VECTORS  TO  WHICH  MVAP  MODEL r 
TO  BE  FITTER. 


BEFORE  SUBROUTINE  RHJONS  IS  CALLED  FOR  THE  FIRST  TIM"  TE 
OBSERVATION  VECTORS  ARE  PLACED  IN  ARRAY  X  DIMENSIONED  NT  BY  II. 
“HE  NC  VARIABLES  TC  BE  PREDICTED  ARE  DENOTED  3Y  FIRST  SUBR  SC'  irT 
•  THROUGH  NC  IN  ARRAY  X  WHILE  SUBSCRIPTS  HOI  THROUGH  NP  DENOTE 
THE  NP-NC  VARIABLES  TO  PF  USED  TO  AID  IN  THE  PREDICTION.  OF  COUP.' 
NP  MAY  EQUAL  NC. 


[*****#*♦«* 


CALL  RH  JONS  (Y,  S,  A„3  1 ,11, S*,  "’C,  LG,  X  M  N,  GAM,  SB,  EE,  DD,  A  B,  B,  BP  ,Q  ,  IT  J 
AFTER  THE  FIRST  CALL  TO  PHJONF,  ARRAY  Si  IS  SCANNED  TD  FIND  THE 
THE  MINIMUM  VALUE  OF  T;!E  AKA  IMF  FPE  PARAMETER  BY  SUBROUTINE 
Hi'-.MIN.  Ty  E  ORDER  MODEL  FOR  WHICH  THE  MINIMUM  OCCURS  IS  LMIN-1. 

CALL  PSC’MTN  (31,LG,LMIN; 

DD  1 1,  ')  =  H 9 9  V S  , 

LG  =  LMI N 

SUBROUTINE  RUJONS  IS  CALLED  AGAIN  TO  DETERMINE  THE  COVARIANCE  AND 
COEFFICIENT  MATRICES  OF  THE  MV  AF  MODEL  WITH  THE  MINIMUM  VALUE 
■  F  THE  AKA  IK  E  FPE  PARAMETER.  DD  (1,  1)  =99999.  PREVENTS  THE 
•^COMPUTATION  OF  THE  LAG-SUM  MATRICES  GAM. 

CALL  PH  JONS  (X,S,  A, SI  ,11,  NP,  NC,  LG,  X  H  N,  GAM,  SB,  EE,  DD,  A  B,  B,  3B,  Q,  ITDT) 
SUBROUTINE  ERRVAR  IS  CALLED  TO  COMPUTE  THE  PREDICTION  ERROR 
MATRICES.  UPON  ENTRY  ARRAY  S  CONTAINS  THE  ONE-STEP  PREDICTION 
E  PROP  COVARIANCE  MAT  Ri  X  AND  ARRAY  A  CONTAINS  THE  COEFFICIENT 
MATRICES.  THE  PREDICTION  F  R  R  OR  COVARIANCE  MATRICES  APE 
■-  TURNED  I  N  ARRAY  BB. 

CALL  EP.RVAP  (A,  B,  S,  BP  ,  NP  ,  LG ) 

STOP 

END 


SUBROUTINE  RHJONS  (X,S,  A,S  1,11,  NP,  NC,LG,  XMN,  G AH,  SB,  EE,  DD,AB,B,BB,2 
1  ,  ITOT) 

V  3  JUROUTI  NE  EH  JONS  COHPDTES  THE  COEFFICIENT  HATBICES  AND  DETER  HINES  IH 
C  PROPER  ORDER  FOB  A  HULTIVARI ATE  AUTOREGRESSIVE  (H VAR)  MODEL.  ON  THE 
C  FIRST  CALL  TO  RHJONS,  NP  TIME  SERIES  EACH  OF  LENGTH  II  ARE  INPUT  INIC 
C  ABBAY  X.  THE  FIRST  NC  TIHE  SERIES  ARE  THOSE  TO  BE  PREDICTED.  LG- 1  IS 
C  THE  MAXIMUM  ORDER  M VAR  PROCESS  TO  BE  FITTED  TO  THE  DATA.  THE  MEAN  OF 
C  EACH  TIHE  SERIES  IS  COMPUTED  AND  STORED  IN  XMN.  THE  AKAIKE  FPE  ARE 
C  COMPUTED  AND  STORED  IN  SI.  AFTER  THE  FIRST  CALL,  SI  IS  SEARCHED  FOR 
C  ITS  MINIMUM  VALUE  AND  THE  INDEX  OF  THAT  VALUE.  ON  THE  SECOND  CALL  TO 
C  RHJONS,  DD  ( 1  (/  1 )  IS  SET  TO  99999.,  AND  LG  IS  SET  TO  THE  INDEX  OF  THE 
C  MINIMUM  VALUE  IN  SI.  AFTER  THE  SECOND  CALL  THE  COEFFICIENT  MATRICES 
C  FOR  THE  LG-1  ORDER  M VAR  PROCESS  ARE  A{NP,NP,2),  A(NP,NP,3),  .... 

C  A{NP,NP,LS).  THE  ONE  STEP  PREDICTION  COVARIANCE  MATRIX  IS  S. 

DIMENSION  GAM  (NP,NP,LG)  ,A  (NP,NP,LG)  ,AB(NP,NP,  IG)  ,  B(NP,NP,LG) 
DIMENSION  BB  (NP,NP,LG) 

DIMENSION  X  (NP,II)  ,XMN  (NP)  ,S  (NP,NP)  ,SB(NP,NP)  ,EE(NP,NP)  ,  DD  (NP,NP) 
DIMENSION  Q  (NC,NC)  ,S1  (LG)  ,  ITOT  (NP,NP) 

DIMENSION  SORK  (  1  000) 

XII=II 

C  H  HEN  DD  (1,1)  EQUALS  99999.  THE  LAG-SUM  MATRICES  GAM  NEED  NOT 
C  BE  COMPUTED.  MISSING  OP.  BAD  DATA  IN  ARRAY  X  IS  DENOTED 
C  BY  THE  VALUE  OF  -100. 

IF  (DD  (1,  1)  .  EQ.  99999.)  GO  TO  123 
DO  1900  1=1 ,  NP 
XMN  (I)  =0. 

ITOT  (1,1)  =0 
DO  1900  J=1,II 

IF  (X  (I,  J)  .  EQ.-  100.)  GO  TO  1900 
ITOT  (I,  1 )  =  I  TOT  (1,1)  *1 
XMN  (I)  =  XMN  (I)  *X  (I, J) 

1900  CONTINUE 

DO  1901  1=1, NP 

XMN  U)  =  XfcN  (I)  /ITOT  (I,  1) 

DO  1901  J=1  ,11 

IF  (X  (I,  J)  .  EQ .  -  100  . )  GOTO  1901 
X  (I,J)  =X  (I,  J)  -XMN  (I) 

1901  CONTINUE 

DO  100  1=1, NP 
DO  100  J=1 , NP 
DO  100  K= 1 , LG 
100  GAM  (I,  J,K)  =0. 

DO  61  L=1 , LG 
DO  61  1=1, NP 
DO  61  J=1  ,  NP 
ITOT  (I,  J)  =0 
DO  62  K=L, II 
K 1  =K-L*  1 

IF  (X  (X,  K)  .  EQ.“100..0R.X(J,K1)  .  EQ.  -  1 00  .)  GO  TO  62 

ITOT  (I,  J)  =  ITOT  (I ,  J)  *1 

GAM  (I,J,L)=GAH(I,J,L)  *X  (I,K)*X  (J,K1) 

62  CONTINUE 

GAM  (I , J, L)  =GAH  (I,  J,L)  *  (II -L*  1 ) /ITOT  (I,  J) 

61  CONTINUE 
123  CONTINUE 

DO  17  1=1, NP 
DO  17  J  =  1,  N  P 
A  (I  ,  J,  1)  =0. 

B (I, J, 1 ) =0. 
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M 


A  8  { I ,  vj  ,  1 )  --  0  . 

BB  {1 ,  J ,  1)  -  0 « 

'?  CONTINUE 

PC  7  1  -  1 ,  N  P 

B{I.I,1)*1. 

A(I,I,1)=1. 

BP  (I,  I,  1)  =1. 

AB  <I,I,1)  =  1. 

DO  7  J=  1 ,  N  P 
S  (I,J)  =G  AH  (I,  J  ,  1) 

SB  (I,  J)  =S  (I,  J) 

7  COSTING E 

IF  (  LG  .LE.  1  )  GO  TC  124 
DO  8  L~  2,  LG 
NL=L- 1 

IF  { NC ,  p. r .  HP)  r;0  TO  2C 
DO  1 1  1-1,  NO 
r'0  21  J-1  ,  NC 

CALL  HATINV  Q,  NC,DET) 

SI  {HI.)  - DET 
20  CONTINUE 

CALL  HATINV  (S,NP,DET) 

IF(NC.NE.NP)  GO  TO  22 
SI  (NL)  =DET 
22  CONTINUE 

CALL  HATINV  {SB,  NT,  DET) 

DO  9  1  =  1, M P 
DO  9  J-  1,  Nu 
DD  (I  ,  J)  =0. 

EE  (I,.1)  =0. 

DO  9  K=1,NL 
K  1  =  L-K ♦  1 
DO  9  11=1, NF 

EE  (I.  J)  =  EE  {I,  J)  -  93  (I,  I1  ,"■>  *GAM  (J,H,K  1) 

DD  (I,. I*  =  DD  (I,  O'  -2  (I,  I  1 ,  E ;  ’’  GAH  ill  ,J,K1) 

‘4  CO  V'T  N  rTB 

DO  VI  7  =  1  ,  NF 
DO  11  J=  1,  NP 
A  (T,J,!.)--0. 

A8  {I,  .7,  L)  =0. 

DO  P  '<  =  1,NP 

A  (I,.J,L)  =A  (I,  J,L)  *DD  (I , K )  *SB  (K,J) 

11  AB  (I.  J,  L)  =AB(I,  J,L)  *EE(I,K)  *S(K,J) 

IF  (L.EQ.2)  GO  TO  12 

DO  1J  K  =  2,  HL 
KN  =  L-K*  1 
DO  13  1  =  1 , N  P 
DO  5  3  J  =  1 ,  N  P 
A  (I,J,K)=B  (I , J , K) 

AB (T, J , K )  =  BB(I,J,K) 

DO  13  K  1=1, HP 

A  (I,J,K)=A  (I,  J,  K)  ♦A(I,K1,L)*BB{K1,J,KN) 

13  AB(I,J,K)  =  AB(I,J,K)  ♦AB(I,K1,L)  *B(K1,J,KH) 

12  CONTINUE 

DO  14  1=1, HP 
DO  14  J=  1,  HP 
DO  14  K= 1 , L 
B(I,J,K)=A (I,J,K) 

14  BB  (I,  J,  K)  =AB(I,  J,  K) 
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DO  15  1=1 ,NP 
DO  15  J  =  1,NP 
S  (I.J)  =GAM  (I,  J,1) 

SB  (I.J)  =  GAM  (I,J,  1) 

DO  15  K  =  2,  L 
DO  15  K1=1,NP 

S  (I, J) =S  (I.J) +A  (I.K1.K) *GAH  (J, K1,K) 

15  SB  (I.J)  =SB  (I.J)  *AB(I,K  1,K)  *GAH  (K1,J,K) 
8  CONTINUE 
124  CONTINUE 

DO  778  1=1 , NC 
DO  778  J=1,NC 
778  £  (I,  J)  =  S  (I,  J) 

CALL  MATINV  (2.NC.DET) 

SI  (LG)  =  DET 
WRITE  (6,1  800)  SI 
DO  200  1=1 , LG 
P1=II^(I-1)  *NP  ♦  1 
F2=II- (1-1 ) *NP-1 
PI  2= (F1/F2) **NC 

200  SI  (I)  =S  1(1)  *F  1 2 
SF  =  31  (1  ) 

DO  201  1=1, LG 

201  SI  (I)  =S  1  (I)  /SF 
WRITE  (6,1800)  SI 

18C0  FORMAT (2X, ’ Si' , 10E12.4) 

DO  77  1=1, LG 

IF  (SI  (I)  .LE.0.)  SI  (I)  =  1. 

77  SI  (I)  =  ALOG  10(S1(I)  ) 

FN=II- (LG-1 ) *  NP- 1 
DO  78  1=1, NP 

DO  78  J= 1 , NP 

78  S  (I.J)  =S  (I.J)  /FN 
RETURN 

END 
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^  ppkmin (si,lg#lminj 

.  '  :  '  J  'E/J  J,i.NE  SEARCHES  THE  VALUES  OF  THE  AKAIKE  FPE  PAR AMET 
JO STAINED  IS  ARHAY  SI.  THE  MINIMUM  IS  FOUND  AND  IHE  ARRAY 
J.  IT-CS  I-  Li.: DIED  bY  LKIN.  THIS  CORRESPONDS  TO  AN  MVAR  MO? 
ORDER  0/  -HI!:-  1. 

DIMENSION  SI  (LG) 

SMIN= 1. E50 
DO  1  I-  1 , LG 

IE  iF.”  (I)  .GT.3HIN)  GO  TO  1 

SHIN  SI  (I) 

COiSJi.i"  .. 

r n  zr:\,z  *(LBjm 

5  FORMAT i2 A, *  LMi  S, FPEMIN* ,  I  i,E13. 5) 

6  FORMAT  riiX,5E*.3. s) 
sett;.* 

END 


SUBROUTINE  HAT  I NV  ( A ,  M, DET) 

THIS  DLUROUTINI  COMPOTES  THE  INVERSE  MATRIX  OF  M  BY  M  MATRIX 
A  AND  RETURNS  THE  INVERSE  IN  A.  THE  DETERMINANT  CF  A  IS 
RETURNED  IN  DET. 

— j.  H  —  «•  s  -  v  .*  h  f  M  ) 

DET"  1 .  C 
DO  1 

PVT  .V  V:,G) 

DET--D-T  “PVT 

a  (  j  ,  j  - ; . 

DO  2  L--  I, H 

.1  a  {J,r,  -a  u.  k’  /vvt 

j  S  1  “  -  *  j  ri 

r-  ■  <r  •  T  1  *1 

«■'  J  /  -*  t  t  f  - 

3  T=A(K,J) 

»i  J)  -  0.  •./ 

DC  «  L-'.,M 

■4  A  {JK,  Li  -A  (K  ,  L)  -  {  A  ( J  ,  L)  *T) 

1  CONTINUE 
RETURN 
END 
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?:  GSR  AM  IS  USED  IC  MAKE  3,  6,  *,  AND  12  HOUR  MVAE  FORECASTS 
D  TUT.'.  TuMPUTE  THE  80*  CONFIDENCE  INTERVALS  TO  EE  PLACED  ABOUT 

am  cghpafe  the  forecasts  i» it.»  the  actual  observations  vali: 

T.i  ?  F  )r  EC  AST  TINE. 

DIN  ION  EMN  (7)  , A  (7,7,10)  ,  A  ED  (7,7)  ,  YDAT  (7,13)  ,  DF.D  (7)  , 

1  I  I  AT  (7,4)  ,XDAT  (7, 9)  ,V  (7,7)  ,  VV  (7,7,  1C)  ,b  (7,7,  1C)  , 

EH  I  (7,4)  ,  FLO  (7,4)  ,  XTQMN  (7)  ,XH:M  N  (7,  PJ  ,  HMCMN  (7,  12) 

D  IT  I.  MIG  N  IVAFN  (28)  ,ISTAN  (21) 

DATA  I S  T  AN/4  H  H  A  N  N  ,  4  HOV  EE  ,  4ii  ,  UHJREM,  4HEN  ,  4Ii  , 

1  4HMCT  2, 4  HINDU,  4H EG  ,4H  PR  AU  ,  4i!NSC;i,  4 ii  WEIG,  4  HK  AOD, 

-i  4 HEOU R , 4HG  ,  4  H  W  E  R  N  ,  4  H  IG  EF  ,  4  UC  u  E  ,  4  .1  WEIS,  4  BEEN  ,4H  / 

FATA  IVAN  N  /4H1ST  ,  4HCLD  ,  4IJLA  YE,  4  UR  iiT  ,  4HC2IL  ,  4HING  , 


1  4  HUT 

,  4H 

,  4HTFMiP,4HEEAT  ,  4iiUR£ 

,4.! 

,  4  H  ,  4  H 

9 

i.  4 

,  411 

,4H/ISI,4HBILI, 4HTY 

,4.1 

,  4R.U-WI,  4 H N D 

9 

3  4H 

,4H 

,4HV-WI,4HND  , 4H 

,  4H 

/ 

KVAE.  Ic  THE  VARIABLE  NUN  EER ...  1  =  H  EIGHT  OF  FIRST  CLC’JD  LAYER  , 
EXCELLING  HEIGHT,  3  =  TENFERATU?E,  5=  VISIBILITY,  6  =  U-hIN'D,  and 
7=  V- WIND.  NT  IS  THE  NUMBER  OF  PREDICT  Zi<  VARIABLES.  LG  IS  TEE 

ORDER  OF  THE  M VAR  MODEL  PLUS  ONE.  IHJS,  FOR  A  91 H  ORDER  MODEL 


LG  WOULD  E  0  U  A  L  10. 

i?IAD(5,2C1)  KVAP  ,  NP,  IG 
iCI  FORMAT  ( 101?) 

L  G  M  1  -  L  G  -  1 
LGMMG-2 
LGF3  =  LG  >3 

THE  MEAN  VECTOR,  CNE-5TEP  PREDICTION  ERROR  CCV 
AND  "HE  COEFFICIENT  MATRICES  FOR  IdE  M VAR  NODE 
HERE.  THE  SAMPLE  MEAN  VECTOR,  WHICH  ID  DE  TERM 
dZFORZ  TEE  LAG-SUM  MATRICES  ARE  COMPUTED,  IS  r 
EMU.  THE  COVARIANCE  MATRIX  IS  BRAE  INTO  ARRAY 


a  COVARIANCE  MATRIX, 
MODEL  ARE  READ  IN 
E  TERMIN'  ED  AND  R  EMCVr D 
IS  FEAT  INTO  ARRAY 
ARRAY  V  AND  THE 


COEFFICIENT  MATRICES  ARE  READ  INTO  Air.  AY  ARD  AND  THEN  PLACED 
IN  ARRAY  A.  A  (I,J, 1)  WILL  ALWAYS  31  THE  IDENTITY  MATRIX  WHILE 
COEFFICIENT  MATRICES  A 1 , A2, A3,...  WILL  BE  PLACED  IN  A  (1,0,2), 

A  { I , J , 3 )  , A  (1,0,4)  , . . ..RESPECTIVELY. 

READ  (5,200)  XMN 
;CC  FORMAT  (7F1  1 .5) 

HZ  AD  (5,  ZOO)  ({V(I .  J)  ,J=1,7)  ,1=1,7) 

DO  5  L= 1 , LG 

READ  (5,  20  0)  (  (AFD  (I,J)  ,J=1,7)  ,1=1,7) 

DC  5  1=1, NP 
DO  5  J  =  1 ,  N  P 
5  A ( I , J , L)  = A  R  E  (I,J) 

IF  Tut  VARIABLE  IS  TEMPERATURE,  THE  GRAND  MEANS,  HOURLY  MEANS 
AND  MONTHLY  MEANS  WHICH  WERE  REMOVED  BEFORE  THE  H V A3  MODEL 
.AS  DETERMINED,  MUST  BE  READ  IN  HERE  SC  Til  AT  THEY  M.AY  BE 
ADDED  cACK  TO  THE  FCFECAST  VALUES. 

IF  (KVAR  .  NE.  3)  GO  TO  110 
READ  (5, 101}  XTCMN 

READ  (5,  10 1)  (  (XHPMN  (I,J)  ,1=1,7) ,J=1 ,d) 

HEAD  (5,  101)  (  (XMCMN  (T,  J)  ,  1=1,7), J=1, 12) 

Id  FORMAT  (7F8 . 3) 

110  CONTINUE 

CALL  LFRV AP  (A, E, V, VV, NP,LG) 

TUI  OBSERVATIONS  TO  BE  USED  TO  MAKE  AND  VERIFY  THE  FOF  ECASTS 
ARE  REAL  IN  HERE.  IT  IS  ASSUMED  THAT  THE  CONTRIBUTION  OF 
THE  HOURLY  AND  MONTHLY  MEANS  HAVE  REIN  nE.IOVE D  BEFORE  THE 
TEMPERATURE  DATA  HAS  3EEN  READ  IN  AND  THAT  TRANSFORMED 
VISIBILITY  AND  CLOUD  HEIGHT  VARIABLES  ARE  TO  BE  READ  IN. 

AH? AY  YDAT  IS  DIMENSIONED  NP  EY  LG+3.  THE  OBSERVATION  TO 
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•  vi  .v-ru  w :::i  thk  12-houp  fofeca.gt  is  ccntainle  ik  llem: 
.  t:,1;,  observation  valid  a:  the  zzec  forecast  tike 

C  .'NT  A  IN  ED  IN  YD  AT  (1,5),  AND  Till  CcS  EE  V'AI  ION  F  AETHZ5T  IN 
.'ADI  15  CONTAINED  IN  YDAT  (I  ,  LG+ 3)  . 

M.A.  (E,  20  ))  (  (YDAT  ( I  ,  J )  ,  I  =  1  ,  N  F  )  ,  J  =  1  ,LGI3) 

J5  1C  1-1, NF 


/ a i  j)  =  yl.m  c,  j*uj  -xxs  o 

Da  2T  I  FOOT  =  1 , 4 

If  (I  FDD  I.  L-l.  1)  GO  TO  21 

L v -  »  —  1  ,  Ni 

4 '  5  .  3  5  -  1 ,  LG  2 

ZDAT  U  ,  LG-  J)  •  X  DAT  (1,1  GX.1-J) 

T  ('.  1;  =FL,V  u  ,IFCST-1) 


..  ...nt 

-  -  ,  •»  .  *•  »i  r 

•  ■  v  4  J 

-  DAT  (I,  IFZi Z)  -  •  • 

DO  4  5  L  -  1 ,  LT1 1 

..  0  4  D  J  -  1  ,  M 

FDA'  (I,  IF  1ST)  -  FLAT  ( I  ,  I  FC5T)  -  A  (I,J,L*1)  ^ M DAT  (J  ,1} 

C  1  NT  I  N” 

1 r  ( X V  A .  .  N  Z .  3 )  GC  TO  21C 

ho -jr  an:  month  cf  the  fop-cas:  zsr.c  tike  ail  i-lad  in 

:  DMY  r  0  N.  T  E.v!  EE  F.  AT ’JR  T  50  THAT  THE  MONTHLY  AND  HO  J  F  L  Y 
S  5  T  AN  3E  ?  D? :  D  BACK  IN.  THE  F  =  1  FOE  DOZ,  InEE =2  FOE  03; 
,i-;ef-.3  foe  did.  ikcf*i  fcf  janjafy  and  tnof  =  u  foe 


-  DA.  (  5 ,  2  5  1)  T  h  ?.  F ,  IXOF 
4 1  1  I--1,N? 

-0  511  J  -  1  ,  4 
I  i:  .3  =  I:i  E  F  ♦  J 

IF  (I.:E.  J7.  a)  IH3=IHP-S 
1  X  C  ~  "  '1  c 

.  ./ A  .  (I,  J)  --FLAT  (I,  J)  +XHFMN  (T,  THE)  ♦MHCMN  (I,IKO)  -XIOMN  (I) 

. i. a .  ; -yd  at  (I,  5-j)  ♦  k  n  (i, the)  ♦  2.1  cuk  (:, : kg)  -xtoks 


^  5  .  I  -  1  ,  N  I- 

3 :  5=1,  u 

FLAT  (I,  J)  -FLAT  (T,  J)  ♦  XN.V  (I) 

-  ‘  J  «.  •  •  »  {i  i  ( ..  ,  4  ,  J )  ) 

:  El  (1,5) M  DAT  (I, J) ♦ 1 . 28* SD 
r 15  (1,5)  -  "DAT  ::,J}-1.28*SD 
15  ML  KV  A? 

:  -  (I„  .GT.  2,  AMD.  ICHK.NE.  5)  GO  TO  3  0 

..  ■  l  :  - 1  c  c  o . 

LF  (T5NK  .  I ..  .  G)  XMUL"=  2000. 

IF  (  r  ■'  T  (I  ,  .1)  .  LE.O  .)  FDAT  (I,J)  =  1  .£-1  5 
I  F  (?:•;:  (I,  J)  .  LD.C  .  )  Fill  (I,  J)  =  1.  E-  1  J 
IF  (1,5)  .LE.O.)  FLO  (I,J)  =1.  r- 1  3 

MAT  C,  J)  =-X.’1ULT*ALOG  (FDAT  (1,0)  ) 

>:  L  C  =  F  LC  (1,0) 

X  ■;  1  =  r  H.  I  C,J) 

F,-;T  (I,  J)  =  -  X  MJLT  *  A  LOG  (XLO) 

"r L ^  (1,0)  =  - XKULT* ALOG  (XHI) 

IF  (FLO  (I, J)  .1/1. G.)  FLO  (I, J) =  C . 
i  DM  (4,  5-J)  =-XMULT*ALCG {YDAT  (I,5-J)  ) 

COM  INN  E 

KV  1=  (KV  AF  -  1)  *4»  1 
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K  v/  -  r*  V  i  ♦  3 

DIE  (0,  9  3  0)  {IVAEN  (KV)  ,KV  =  KV1  ,NV2) 

}C  F'Jr.UAI  (/2X  , 4A4/) 
do  *♦  ^  i  =  i ,  n  r 
isn=  i:-i)  *3*i 

I^r^=IST1 +2 

WPITL  (t> ,  40  1)  (1ST  AN  (1ST)  ,  IST=  I  ST  1 ,  1  3T  2) 

:  1  FOF/iA".'  (/2>:,  3A4/) 

*717;.  (c,4l)  I  ,  (FLAT  ( I ,  J)  #  J= 1 #  4 )  /  i  l  0  A  .  ( <■  t  j  /  >-  -  1 »  L‘ ) 

41  FORi-AI  Uv,  '  !=•  .15,'  FCST* ,4El 3. 5, *  ACT  0  A  L  '  ,  E  1  3  .  -> ) 
w:i;e  (d,42)  (flc  (i,J) ,j=i ,4) ,  (fh:  (i  ,J) , j=i , 4) 

42  ruSSAT  •  .S3  LC*  ,  4E  1  3.  5,  •  .S3  -i  I  *  ,4213.3; 

4  3  C  0  4  7  7'.  'J  E 


2  OF  DOLT  IN  2  EFF.VAR  (A,B,V,  V7,N?,L3) 

nE.iii.sio:;  V{NP,hp)  ,a  (NF,  np, lh)  ,n  (np,  np,  lg)  ,  v v  (nf ,np, LG) 
DC  1  1  =  1,  N'P 
DO  1  J= 1 , N  P 
1  VV  C,J,  1)  -  v  (I, j) 

LG  E  1  =  LG  -  1 

DC  r,  L=1,LGM 

DC  o  2= 1 , N P 

DO  •)  J  -  1  ,  !J  P 

DO  ^  L  "■'  =  1 ,  L 

;  ?  (Lh.GT.  1  )  GO  ?C  7 

t  ( 1 ,  J  ,  L )  =-A  (I,J,L*1) 

gc  :  i,  !, 

7  ccnt  t 

8  J \:  ,j’/L)  =  y  (I,J,L)  -E(I,KK,LK-1)  *A(KK,J,LK) 
u  CONTINUE 

DO  10  L=  2 , LG 
DC  13  I  =  1  ,  N  P 
DO  1C  J  =  1 ,  N F 
VV  (I,J,L)  *VV  (I.J.L-1) 

DO  10  E  K  =  1  ,  N  ? 

DO  10  LL=1,KP  14 

r  VV  (1,0,  L)  =  VV  (I,J,L)  *9  (I,KK,L-1)  *VV  (Erl,--,  1)  *B  1) 


FLOWCHART  OF  ANALYSIS  PROCEDURE 


Initial ization 

f - ! 

|  Before  any  subroutine  ceils:  II  is  the  number  of  NP-  j 
I  dimensiona1  observation  vectors  an  MVAR  model  is  to  ! 
i  be  found  for,  the  maximum  order  model  to  be  tried  is  I 
|  lG-1,  the  array  X  dimensioned  NP  by  II  contains  the 
;  observation  vectors.  I 


1st  Call  o *  Subroutine  RHJONS 


vvAR  "iode":  >  Tr  order  zero  to  order  LG  are  fitted  to 
tne  II  N°'--u  imensiooal  observation  vectors.  The 
values  of  too  Akaikc  r':E  carameter  are  stored  ’n 
array  Si. 


Call  Subroutine  FPEMIN 


■  Array  SI  is  searched  for  the  minimum  value  of  the 
I  FDE  parameter.  Array  element  LMIN  denotes  the 
|  minimum  ana  LMIN- 1  is  the  order  of  the  MVAR  model. 


2nd  Call  of  Subroutine  RHJONS 


An  MVAR  o-ode’  of  order  LMIN-1  is  fitted  to  the  II 
j  NP-dimersional  observation  vectors.  The  one-step 
prediction  error  covariance  matrix  is  contained  in 
|  array  S  and  the  coefficient  matrices  in  array  A. 


Cali  Subroutine  ERRVAR 


j  nredict'on  error  covariance  matrices  for  various 
j  stec  forecasts  are  computed  and  contained  in  array 
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